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ANALYSIS SYMBOLS* 


exponent for a basic power-law camber line, see Eq. 36a 
see Eqs. B-10 


Fourier coefficients for the trailing jet-sheet vorticity 
distribution for a singularly blown flat plate, see Eqs. B-9 


wing span 
see Eqs. B-10 


Fourier coefficients for the trailing jet-sheet vorticity 
distribution for a regularly blown flat plate, see Ref. 4 


airfoil chord (taken as c = 1) 
flap chord 

jet-momentum coefficient, fis/a,.¢ 
lift coefficient, (0) + oft) 


(0) 


interference lift coefficient, Cy - c,° = he, 


lift coefficient without jet-flap blowing, i.e., for unblown 
airfoil 


lift-coefficient partial derivative, de, /aa 
lift-coefficient partial derivative, de, /38 
lift-coefficient partial derivative, dey [dr 
lift-coefficient partial derivative, de, /dK 
pitching-moment coefficient about the leading edge, positive 
t96 a bees 2 9 aa moment (also see alternative definition below), 
of) + c 
quantity defined by Eq. B-11 (also see alternative definition 
above ) 

) 


interference pitching moment coefficient, o - 0° = dcy 


“Computer program symbols are listed in Appendix G. 
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ANALYSIS SYMBOLS (Contd) 


moment coefficient without jet-flap blowing, i.e., for unblewn 
airfoil 


moment-coefficient partial derivative, dc,,/ 30. 
moment-coefficient partial derivative, ac / 06 
moment-coefficient partial derivative, dc / chs 


pressure coefficient, (p-p_)/q,3 for small perturbations, 


cy = -2u'/U,, 


Fourier coefficients for the trailing jet-sheet vorticity 
distribution for a power-law cambered jet-flapped airfoil 
at zero angle of attack and zero jet deflection (+ = 0), 

see Eqs. 49 ard 52 


see Eq. B-11 


Fourier coefficients for the trailing jet-sheet vorticity 
distribution for a regularly blown mechanically flapped 
airfoil at zero angle of attack, see Eqs. B-15 and B-19 


influence function coefficient in quadrature method, 
3°D,/ae?, see Eqs. C-1 and C-3 
defined by the generalized set of simultaneous equations, 


Ent dan) En = ent fq» where the coefficients E,, e,, and f, 
n= 


may be those corresponding to any one of the boundary-value 
problems corsidered in this report, i.e., E, may be A, Bye Cy» 
» E,, ete. along with the appropriate analytical relation for 


D, 
en and f. (used principally in the computer program) 


Fourier coefficients for quadrature-method influence 
coefficients, see Eq. C-4+ 

see e. above and tn be low 

see Eq. B-20 


see Eq. 51 
see Eq. 42 


nondimensional vorticity distribution along the trailing jet 
sheet, ¢'(x) = ¥(x)/U,¢ 
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see Eq. 50 

maximum camber height, see Fig. 2 
see Eqs. 13 and 28 

see Eqs. 28 and 30 


see Eqs. 28 and 31 
see Eqs. 28 and 32 
see Eqs. B-4 

see Eqs. B-6 and B-7 
see Eqs. 26 and 29 
see Eqs. 29 and 33 
see Eqs. 29 and 34 
see Eqs. 29 and 35 
see Eq. C-2 

see Eq. D-12 

see Eqs. 44 

see Eqs. 3-8 

see Eqs. D-10, D-14, D-22, and D-26 


meh collocation point, m = 0,1,2,. _ .,N-1 


see Eqs. D-11, D-14, D-23, and D-26 
momentum flux per unit span within the jet sheet 
degree of polynomial for a polynomial camber line, see Eq. 37 


total number of collocation points {also total number of 
simultaneous equations) 


summation index, n = 0,1,2,. . - 
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ANALYSTS SYMBOLS (Contd) 


P static pressure 
a, freestream dynamic pressure, (p/2)U2 
r see Eqs. 36c 
s see Eq. C-6 
S wing area 
s A see Eq. 23 
Sp see Eq. 24 t 
Sh influence coefficient for quadrature method, 3°s_/as’, see Eq. 27 
T aircraft thrust | 
U,V local velocity components, respectively parallel and ; 
perpendicular to the freestream velocity + 
u,v local perturbation velocity components, u' =u - Us v'sev 
U. freestream velocity | | 
U (x) unit step function;€@((x) =0 for x < 0,&/(x)=1 fox x >0 
j W aircraft weight 4 
Xs rectangular coordinates, see Figs. 2 and 3 
y airfoil camber line ordinate 
[1 = (tex) BY/L1 4 (1=x)9) | 
a angle of attack | 
Oe, angle of attack at zero Lift | | 
Bh "equivalent" Fourier coefficients for quadrature method, 
see Eqs. 10 and 11 | 
Y vorticity | 
i, ay/aa 
%, ay/36 | | 
| ’. ar/ar 
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ANALYSIS SYMBOLS (Contd) 


flap deflection angle, see Fig. lc 


interference lift coefficient, c : - of?) = oft) 


interference pitching-moment coefficient, cy oe off) 
A, 7 Dy» see Eq. 12 

Sa7Sps see Eq. 22 

dummy variable, see Eq. 38 

local camber-line slope (also see below) 


polar coordinate defined by x=(#)(1+cos6) for O0<x<1, 
see Eq. B-16 


camber ratio, h/c 

4/e j 

dummy variable along abscissa x, see Figs. 2 and 3 

location of flap hinge point, £ = (#)(1+4cos%) 

pivotal points for quadratic function approximation in an 
infinitesimal region at the airfoil leading edge; 0 < Eo < & 
pivotal points for quadratic function approximation in an 
infinitesimal region at the airfoil trailing edge; Eo < §3 <1 
and 4 7 1 

fluid density 


jet-sheet deflection angle at the airfoil trailing edge, 
measured relative to the airfoil chord line as depicted in Fig. 1 


polar coordinate defined by x = cos~“(p/2) for x >1, see Eq. 45 
mi /N3 m=0,1,2,- ~ .,N-1 

polar coordinate of flap hinge point, defined by € = (4)(14 cos X) 
denotes the airfoil leading edge 

denotes the airfoil trailing edge 


denotes the trailing streamline (i.e., the jet sheet for the 
blown flow case) 
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SECTION I 


INTRODUCTION 


The most commonly proposed application of the jet~flap principle to 
date has been to jet-augmented mechanical flaps for use during the take-off 
and lending flight phases of aircraft. For this application, the lift 
coefficients achieved are very high and the reletive lift contribution 
due to camber-line jet-sheet interaction is negligible. Recently, studies 
have been conducted on the aerodynamics of pure* jet~flapped wings at high 
subsonic and transonic speeds (Refs. 1 and 2). Motivation for these 
studies is the possible use of pure jet-flapped wings for maneuvering 
combat aircraft at high speeds. For pure jet-flapped wings at high speeds, 
the camber~line jet-sheet interaction, although small in absolute magnitude, 
can be a relatively higher fraction of the total lift than in very-high-lift 
applications and should be taken into account in some instances. Because 
Of previous emphasis on very high-lift applications, little attention has 
been given to the analysis of camber-line effects. Other than a very 
limited investigation by Woolard (Ref. 2), the sole in-depth treatment of 
camber-line effects appears to be that of Hough (Ref. 3), who formulated 
an analysis for a polynomial camber line, but provided specific numerical 
results only for the parabolic shape. There is a need for a prediction 
capability for arbitrary camber-line shapes, since the parabolic camber 


line is not necessarily the best one for high-speed applications. 


The analysis herein is a small-perturbation one with the obvious 
restriction to small flow disturbances and the attendant results that 
flow solutions and boundary conditions are linearly superposable. Hence, 
the flow about an arbitrarily shaped jet-flapped airfoil may be obtained 
by appropriate addition of jet-flapped-airfoil solutions for a regularly 


"That is, a wing employing trailing-edge jet blowing alone, unassisted by 
mechanically deflected flaps or ailerons. 


l 
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blown* flat plate at angle of attack (Fig. la), a singularly blown** flat 
plate at zero angle of attack (Fig. 1b), a singularly blown cambered plate 
at zero angle of attack (Fig. 2), and the solution for the flow about an 
unblown symmetrical airfoil at zero angle of attack. To the order of 
linear terms in small perturbations, airfoil thickness does not influence 
the vorticity distribution nor the associated lift and pitching-moment 
coefficients. The local surface-pressure coefficient is affected, however, 
and is the sum of contributions due to the vorticity and thickness 
distributions. In view of the superposition properties of solutions and 
boundary conditions, the present analysis is concerned principally with 
camber-line aerodynamics on the assumption that the solutions for the 
other aforementioned contributions are available in the technical literature. 
The regularly and singularly blown flat-plate solutions, for example, are 
available in Ref. 4. The singularly blown solution is given also in an 
Appendix in this report. Solution for the flow about a symmetric airfoil 
of arbitrary thickness distribution at zero angle of attack may be found 

in Ref. 5. The camber-line solutions obtained herein are, of course, 
limited to camber-height ratios which are small relative to unity as a 


consequence of the small-perturbation assumption. 


Consideration of high-speed aerodynamics herein is limited to sub- 
eritical flight Mach numbers so that the jet-sheet camber-line flow- 
interaction problem can be solved for incompressible flow and extended 
into the compressible regime by means of jet-flapped airfoil similarity 


rules such as given by Woolard in Ref. 2. 


In this report the Hough formulation for polynomial camber lines is 
further developed (with modifications) into a method designated as the 
"Power-Law Superposition Method." Since it is likely, however, that a 
polynomial representation may be inadequate or unwieldy for some camber 
lines, a new method designated as the "quadrature method" is develope¢ 


* Regular blowing denotes a jet sheet emerging tangent to the trailing edge. 
**Singular blowing denotes a jet sheet emerging at an angle relative to the 
trailing edge. 
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to avoid the use of a polynomial and permit the use of numerically 
specified ordinates for completely arbitrary shapes (except for slope 
limitations at the leading and trailing edges). 


The quadrature method is the principal contribution of this report 
and is the method emphasized. This method is derived from superposition 
of jet-augmented, mechanically flapped airfoil solutions (Ref. 6) taken 
in the limit as the number of flapped airfoils become infinite. The 
approach ultimately yields the aerodynamic properties in terms of integrals 
having integrands which consist of the product of the camber-line ordinate 
and an influence function. In application, the integrals involved are 
evaluated by numerical methods. The required influence functions are 
determined by means of a series-collocation solution to the governing 
integral equation. Tables of the influence functions, obtained from 
nine-point collocation solutions on a CDC6600 computer, are presented in 


the report. 


The Power-Law Superposition Method for an arbitrary camber line 
employs superposition of (n-1) solutions for basic power-law camber-line 
shapes of the form r,(x-x"), where @ = 2,3,~----- >n. This approach differs 
from that of Hough who uses superposition of camber lines of the form 
a The analytical bookkeeping involved in studying a particular camber 
line is believed to be more convenient with the present approach. The 
governing integral equation for each basic power-law camber line is solved 
by collocation assuming a Fourier series approximation. Numerical 
solutions for the Fourier coefficients were obtained on a CDC6600 computer 
for nine collocation points, seven basic camber-line shapes corresponding 
to a = 2,3,4,5,6,7, and 8 and values of C, ranging from 0.001 through 5.0. 
Tables of the coefficients so obtained are presented in the report. 


One of the principal features of the quadrature method is that it 
provides a convenient means for determining the aerodynamics of an arbitrarily 
cambered jet-flapped airfoil by hand calculation utilizing a pocket- or 
desk-electronic calculator. The hand calculation approach is very 


t 
| 
t 
i 
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convenient* for the lift and moment coefficients, but somewhat less so 

for pressure distributions and trailing jet-sheet shapes. If, therefore, 
one requires aerodynamic characteristics for an arbitrary camber line for 
only a moderate number of parametric variations, hand calculation by the 
quadrature method is appropriate. If, on the other hand, a large number 
of parametric variationsis involved, the sheer bulk of the calculation 
dictates the use of a digital computer. For application convenience in 
this latter situation, or for the purpose of calculating additional tables 
of influence functions, a digital computer program is given in an appendix. 
If the airfoil camber line is closely represented by a polynomial, the 
power-law superposition method is likely to be the most convenient. 


For analysis convenience, a given aerodynamic property of a jet-flapped 
airfoil with the jet operating (blown airfoil) is expressed herein as the 
sum of the aerodynamic property without the jet operating (unblown airfoil) 
plus an interference quantity; that is, blown-airfoil aerodynamics equals 
unblown-airfoil aerodynamics plus an interference quantity. The subject 
of this report is the determination of the interference quantities for the 
principal aerodynamic properties. In application it is assumed that the 
detailed methodology for calculating the theoretical aerodynamics of the 
unblown airfoil is obtained from one of the many sources in the technical 
literature (Refs. 5 and 7, e.g.). Nevertheless, for convenience, a brief 
summary of some principal results from unblown thin-airfoil theory is 
presented in Appendix A. 


As mentioned earlier, the most likely application of cambered jet- 
flapped airfoil theory is to pure jet-flapped wings on combat aircraft 
maneuvering at high subsonic speeds. For this application, it is of 
interest to estimate an upper bound to the value of the two-dimensional 
jet-momentum coefficient. Employing a mean wing chord, (S/b), a mean 


two-dimensional jet-momentum coefficient, c,, may be cast in the form 


e" 


c= kg(W/S)(T/W)/q., 


*Especially for the small personal programmable calculators 
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where (W/S) and (T/S) are the wing loading and thrust-to-weight ratio 
respectively and k, is a factor which is proportional to the fraction of 
engine thrust available to the jet flap, the fraction of sea-level thrust 
available in flight at altitude and an engine installation factor. 
Assuming a wing loading of 80 pounds per square foot, a thrust-to-weight 
ratio of unity, a maneuvering Mach number of 0.9, and a maneuvering 
altitude of 25,000 feet, to be typical for a fighter aircraft, and very 


conservatively assuming k 


to be 0.5 yields €; = 0.10. Hence, it appears 


f 
that for foreseeable high-speed applications, values of the jet-momentum 


coefficient that are of practical interest are less than 0.1. However, 
for the sake of completeness, and because future applications are difficult 
to anticipate, the tables and examples presented in this report include 


values of e, ranging from 0.901 through 5.0. 
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SECTION II 


ANALYSIS 


THE QUADRATURE METHOD 


As noted in the introduction, the aerodynamics of a cambered jet- 


flapped airfoil (Fig. 2) by the quadrature method is derived from super- 


position (Fig. 3) of solutions for regularly blown, mechanically flapped 


airfoils taken in the limit as the number of flapped airfoils become 


infinite. In the present analysis, the boundary condition of zero jet 


deflection at the trailing edge (Fig. 2) for the camber case is evoked 


to achieve, for the more general airfoil geometry ana aerodynamic state, 


a superposition rule in which the contribution of each fundamental case 


is proportional only to the single principal parameter for that case. 
For example, for the lift coefficient, this approach yields the relation 


%* %." + 5 Sadad ~" + 2% (1) 


with corresponding relations for the other aerodynamic properties. If 


tangential trailing-edge blowing (regular blowing) were to be employed 


as a boundary condition, the singular-blowing contribution would be 


proportional to (6.5 - T) rather than to T, and the magnitude of Cp 
kK 
would differ from that in Eq. l. 


Before developing the quadrature method, it is of interest to discuss 
a pertinent limiting property of Spence's solution (Ref. 6) for the 
regularly blown, mechanically flapped airfoil. Consider the fundamental 
cases of the regularly and singularly blown flat plate (Ref. 4) and the 
regularly blown, mechanically flapped airfoil illustrated in Fig. 1. In r 
order to properly implement the quadrature method as formulated in this 


report, the aerodynamic derivatives cp (c,+6), cn (e556), and Valo, »x58) 
for the regularly blown, mechanically flapped airfoil should reduce to the 
corresponding derivatives with respect to a for the regularly blown flat 
plate and with respect to T for the singularly blown flat plate for the 


limiting cases of € = 0 and 1 respectively. These limiting cases are 
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examined in Appendix B where it is found that the limit for — = 0 yields 
the proper result, but that for § = 1 does not. Hence, for the purpose 
of distinguishing between quantities resulting from taking € = 1 in the 
Spence theory and those obtained directly, the Spence quantities will be 
denoted by an inverted circumflex, e.g., ep (c, rar a (c, 31), and 
Yale, 2.1). . . 


In order to overcome the aforementioned difficulty at € = 1, the 
vorticity-distribution derivative with respect to 6 must be written as 


My lre€) = Mane) + Ene) 4-1) -H] cay 


where Y is the unit step function for whichM(E-1) = 0 for (E-1) < 0 
ana Y (€-1) = 1 for (E-1) > 0. For the sake of conciseness of presentatio: 
the dependence of Yg and Y, on es is not symbolically indicated. 


From the principle of superposition, the vorticity-distribution 
derivative for the segmented camber line of Fig. 3 is 


N-l 
lx) = - Flo) M 2 % (x, €,)(On-On-1) + lx On-1 (3) 


where 


708, = 12, = = 0, and 


N 
t= 2 (0, -G_,)=0 


ase 
Noting that 
@ = dy/d€ sy’ (4) 
and hence 
(8n-On-1)= AO= Ay’ = (Ay/AE)AE (5) 


where AE ZTE -€ Substitution of Eqs. 4 and 5 in Eq. 3 yields 


n n-1° 


N-I 
Y(x) =— %(x,0)y"(0) ~2 Ty (ts€q)(Av/Q€) OE ty (eIde(H) (6) 


T 


Sine ah Sd 
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where y'(1) is the slope of the airfoil camber line at the trailing edge. 
Taking the limit AE + 0 in the second term on the right-hand side of 
Eq. 6 yields 


rond=[re CeCe] = fel EGAE VE cr 


Substituting Eq. 2 in Eq. 7 yields 
(0) vi) 
(x)= 7°" x) & (i) Men) (x,1)] 


+ [w(é %(n,e)], - fuse) ei/ae* ae (8) 


where, assuming that (ay/aE*) is bounded at the trailing edge, the 
integral involving UY (§-1) vanishes. 


Assuming continuous first and second derivatives for Ys and y and 


integrating Eq. 8 twice by parts yields 


(i) 


(x)= 7a) + (1) [%0)-% (xs1)] 
- fale Na ae?) ag 
Substituting Eqs. B-1 and B-15 in Eq. 9 yields 


Y(1)/Ue = 7 x)/Ue -(2/m)y"(1) 2% fal 1x) 


+[4x/( + x)] PAR + 2x°* < XB, 


ee 
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where X is defined in the list of symbols and 


Bayt ADn- Hp (11) 

ADp (cj) An (c)) - Dple,,1) (12) 
i} 

Hy = ff vl€) 0% (c,,€)4€ (13) 
° 


The determination of AD (c;) is treated in Appendix B, with tabulated 
numerical values given in Table I*. For the sake of completeness, 
tabulated numerical values of A,(e,) are presented in Table II. In 
general, the integral in Eq. 13 must be evaluated by numerical means. 
Discussion of the methods for accomplishing this is deferred to a later 
point in the report. The determination of Dile,»&) is treated in Appendix 
C, with tabulated numerical values given in Table III. 


For zero jet deflection at the trailing edge (t = 0), the lift 
coefficient and the tail-down pitching-moment coefficient about the 


leading edge are given, respectively, by 


| 
\ 
%° 2 [o/-) (14) 
' 
Cem = 2 fx(y/ve dx (15) 


It would appear that the logical procedure for obtaining cp and c. 
in terms of the 8 coefficients would be to substitute Eq. 10 in Eqs. 14 
and 15. This is the only option available for cw but a simpler relation 
can be obtained for Cp by an analysis paralleling that for Eqs. 97 through 
103 in Ref. 4. Noting that Eq. 10 is of a form similar to that of 


*Since Tables I through IV and VII and IX are reproduced from computer 
printouts the entries in the tables are identified by the Fortran symbols 
for the parameters listed. The Fortran symbols and their corresponding 
analysis symbols are listed in Appendix G. 
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Eq. 107* in Ref. 4, then by analogy with Eqs. 103 and 107 in Ref. 4 for 
the lift coefficient and by direct substitution of Eq. 10 in Eq. 15 for 
the moment coefficient, there is obtained 


_ (0) 6 
=o + 4B, (16) 


N-1 
Gm = Om ~¥'NIp- EL Bala (17) 
n=0 


where I, and I, are given by Eqs. B~4, B-5, and B-6. 


It should be noted that in view of the similarity of Eq. 10 to 
Eq. 107 in Ref. 4, and other mathematical similarities (without exact 
equivalence) to Ref. 4, some of the mathematical relations given herein 
May be deduced by means of a term-by-term comparison with those in Ref. 4 


without recourse to actual mathematical manipulation. 
By comparison with Eq. 123 in Ref. 4 (taking into account differences 


in coordinate-system sign convention), the trailing streamline (jet-sheet) 
shape is given by 


(o) .., - 
¥45(?) = a -2y'(!) (1- cos ) -2B,| 4 tan (= + m)-sin =] 


-5 By 2(cos + sinn? -2nsin < cos np) / (an? - 1) (18) 


x= cos “(/2) (x21) 


*When Eq. 107 of Ref. 4 is corrected for significant typographical errors. 
See Appendix B. 
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The derivatives c, (c,,&) andc_ (c,,&) may be modified in a manner 
] es Mm. j 
analogous to Eq. 2 so that 


cag(eve = eg (6) + 4m [Dg(c),€) +4 (€~1)A09(c,.€)} (20) 
em(cj>6) = ems (€)- Sp(¢j€) -M(é-1) [19 + As(c,,€)] (21) 


where 
N-| 
As ® S - S) = YADgIp (20) 
n=O 
N=] -| 
he D Anta Sp = 2 Ontn (23) (2h) 
n=0 n=0 


Tabulated numerical values of As(c,) are given in Table I. 


Making use of Eqs. 20 and 21, relations for the lift and moment 
coefficients also may be derived by a superposition procedure paralleling 
that used to obtain Eq. 10. The resulting lift coefficient is identical 


to Eq. 16. The moment coefficient, however, is of a different form given 


by 
Cm= cin) = yy + As) + SY (25) 
where 


9 = [v(e)s5(c).€¢€ (26) 


$5 (c;.€) =F 056) Tn (27) 


and tabulated numerical values of Sple 578) are given in Table III. 
Determination of the pitching-moment coefficient by means of Eq. 25 


involves less computational labor than does Eq. 17. 


CT Stn me  * 
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The integrals of Eqs. 13 and 26 must be evaluated by numerical 
quadrature. However, as may be seen in Appendix C, the function pie 
and hence also Sp? experiences a singular behavior at — =.0 andl. It 
follows therefore that the leading and trailing edges must be given 
separate special treatment in the aforementioned integrals. For this 
purpose it is convenient to define the following 


(c) 
Hy = (Ha) e + i. + (Hn). 


,) : I, ; & - 


TE 


é, fo 
(Hn), ~ Jol€)on (€)ag ; Ha” = [veepon(edee (30) (32) 


Cr) es J le)os(eyee 32) 


é " é 
De = [ole sale )a€ Ss [lesi(edet (33)(34) 


Se y(€)sn(€)ag (35) 


with & << 1.0 and (1-€,) << 1.0. Approximate analytic relations for the 
foregoing leading and trailing edge integrals are given by Eqs. D-8, D-9, 
D-20, and D-21 in Appendix D. Evaluation of these analytic relations 
requires values of ECE) and s(&) at five selected magnitudes of &, 
which were taken to be 0.025, 0.050, 0.975, 0.9875, and 1.0000 in this 
report. Tables of values of ECE) and s(&) for these specified arguments 
are given in Table III. 
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The integrals of Eqs (31) and (34) are now the ones to be evaluated 
numerically. The specific computational algorithm employed for this 
purpose will depend upon the precision required, the type of computational 
mechanics used, the character of the quadrature integrand, and perhaps, 
the number of cases to be calculated. For purposes of discussion it is 
convenient to categorize the computational mechanics as either a 
"nonautomatic hand computation," an "automatic hand computation," or a 


"machine computation." 


A nonautomatic hand computation employs non- 
programmable calculating machinery and requires that the airfoil camber- 
line coordinates and the influence function be input by hand at intermediate 
stages during hand manipulation of the computation steps on the calculator 
keyboard. An automatic hand computation employs limited-capability 
programmable calculating machinery* and, because of the machine limits, 
requires that the camber-line coordinate and influence function be input 

by hand at intermediate stages during the automatic running of the 
computation steps on the calculator. Both nonautomatic and automatic 

hand computations require tables of the influence functions Die; .8) 

and Site, »6)- A machine computation employs a digital computer, permits 

all camber-line coordinates to be initially input together and calculates 

or stores the influence functions internally to the program. The quadrature 
algorithm complexity and integration step size that can be tolerated 
obviously are strongly influenced by the type of computation mechanics 
employed. 


The quadrature integrands in the present application vary rapidly 
over the aft portion of the airfoil and also are of a greater magnitude 
in this region. In selecting a numerical quadrature algorithm for good 
accuracy, therefore, one has the choice of employing a large number of 
equally spaced abscissae or a lesser number of unequally spaced abscissae 
with the spacing graduated so that the spacing decreases over the aft 
region of the airfoil. A nonuniform abscissae spacing can be achieved 
by dividing the range of integration into subintervals each of which 
employs an equal-abscissa-spacing quadrature formula, but for which the 


*Such as, for example, the Hewlett-Packard models HP-25, HP-67, and HP-97. 
or the Texas Instruments models SR-52 and SR-56. 
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spacing differs for the various subintervals, or by deriving specific 
formulae for nonequally spaced abscissae. The use of a large number of 
equally spaced abscissae increases the labor required to input the data 
regardless of whether the computation is performed by hand or machine. 
For hand computation, however, this approach requires more extensive 
influence-function tables if the tedious process of interpolation is to 
be avoided. Dividing the range of integration into subintervals also 
may introduce difficulty in achieving abscissae values coinciding with 


those listed in the influence-function tables. If only a few integrations 


are to be performed, use of a large number of abscissae, the subinterval 


scheme, or even influence-function table interpolation may not be 


objectionable. 


For nonautomatic hand computations it is likely that, in general, 
the most practical quadrature formula is the trapezoidal rule with 


nonequal abscissa spacings graduated in the manner previously noted. 


Although a computer program for machine computation is included in 
Appendix H, the principal emphasis is upon automatic hand computation. 
This is largely because of the wide availability of programmable pocket 
calculators. The capabilities of the key-programmable Hewlett-Packard 
HP-25 and Texas Instruments SR-56 models were particularly influential 
in establishing the level of complexity of the quadrature formula selected 
and presented in Appendix E. The formula of Appendix E employs a quadratic 
interpolating polynomial for the integrand, arbitrary nonequal abscissae 
spacings and an arbitrary number of abscissae. Program listings and run 
instructions for use of this formula on the HP-25 and SR~56 calculators 
are given in Appendix F. 


Although the selection of the aforementioned quadrature was dictated 
by considerations of use in automatic hand computations, for convenience 
and consistency the same formula is also used in the digitel computer 


program. 
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THE POWER-LAW SUPERPOSITION METHOD 


The "power-law superposition meth employs superposition of basic 


power-law camber lines of the type 
¥y = fq (xx) a= 2,3,4,--- (36a) 
for which the abscissa for the maximum camber height, h, is 


= o ten) (36b) 


and f. is given by 


to = Kq (eye xh) (36c) 


For an arbitrarily shaped camber line given by 


M 
y= L ytx-*) (37) 
a=2 


the Ke coefficients are determined by the specifications placed on the 


camber-line geometry. For example, for an S-shaped camber line of amplitude 


h, and antisymmetric about the mid-chord point, the coefficients are 


ry 1873 h and 2,” -1273 h. 


From Eq. 69 of Ref. 4, the integro-differential equation for the 
vorticity distribution, g'(x), along the trailing streamline is 


£ woe (E(B 


2 sx-1\'42e! Ve a7) 
-+(#) 8 (=) =i 7 = 
where g(x) = (1/2)e s¥ng (x) and the boundary conditions are 


g(t) =g(@)=0 (39) 
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The chordwise vorticity distribution along the airfoil is given by Eq. 62 
in Ref. 4 as 


(x) yloh - o 
te oe a(S) SG) EP to 


In the process of solving Eq. 38 it is convenient to treat separately 
the logarithmic singularity in the vorticity distribution (Ref. 4) known 


to occur at the trailing edge. For this purpose, using the form of Eq. 2.26 


in Ref. 3 as a guide, let 


9(x)= [ 2rq(o-10 m7] F(x) + (172) q(x) (2) 


where 
Fox) = -(2/x"2pen(x-t) +2La[(x% (x24 0] (42) 


with G(1) and F(1) finite and G() and F(~) vanishing. 


Substituting Eqs. 36 and 41 in Eq. 38 yields 


2 p® 
+ (*) La dn -XG(x) 


~ al aj \ 2 
22 {292 Y Fi ~ 247 +aicose| (22 - | 


Pek ee 
-20(—)  Kj,Xx 
( x ) im 4 (43) 
where 
Kgs [1-3-5 (24-1]7 (2-4-6 2a) (44) 
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For the purpose of solving Eq. 43 it is convenient to introduce 


the transformations 
x=cos*(p/2); 7 = cos *(8/2) (45) 


and expand G'(®) in the Fourier series 


eo 
G' (@) =8in $ 2 on cos np (46) 


in the range 0 < @< 7. The corresponding series for G(~) is 


ica 
? ; : 
Gip)= > 2C,, (cos cos np + 2n sin = sin np )/(4n?-1) (47) 


n=O 


Substituting Eqs. 45, 46, and 47 in Eq. 43 and satisfying the resulting 
equation at N points given by 

@m=mm/N, m=O0,1,2,----,N-I (48) 
yields the result 

N-I - 
ry (Omn+ ADmn) Che frm + \Gm /4 (49) 


where an and Dan are given by Eqs. A-8 and 


= [ £a(ton Em) ~ sec Em fn(ton <m)| (50) 


im 28 (0-1) cos® 2m + (See) 
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A word of caution is in order regarding precision in the numerical 
calculation of * by Eq. 51. For a given value of a, as m~*N, small 
differences between small numbers are encountered in the bracketed terms 
in Eq. 51. The situation is greatly aggravated fora~+N. If, therefore, 
a given precision in a is desired, care should be taken to insure that 
the bracketed terms are calculated to a sufficient number of places to 


achieve the desired precision. 


The N simultaneous equations for . given by Eq. 49 have been 
numerically evaluated on a CDC 6600 digital computer for N = 9, a = 2, 
3, 4, 5, 6, 7, and 8 and values of c, ranging from 0.001 through 5.0. 
Values of th + rel were input to 10 significant figures* and the 
calculations performed to 15 decimal places. The resulting values of 
c are recorded in Table IV. The computer program for accomplishing the 


preceding computations is documented in the listing in Appendix H. 


Making use of Eq. 47 in Eq. 41 and substituting into Eq. 40 the 
value of g'(n) obtained by differentiating Eq. 41 yields the following 
for the vorticity distribution along the airfoil 


Y (4) /U @ = y nu, +rgla-I)2/m)x *“gn(1-x) 


N-I 
+ (19/4) [4cgX+X) +25 Cy x"| x 3/2 
n=l 


where X is defined in the list of symbols. 
By analogy with Ref. 4, there is obtained 


cg= cg + rqCo 


N-I 
Gm = Cm'= ral (1-0) 14 + 17a Sata] 
= 


*Because of the previously noted precision problem, double-precision 
arithmetic was employed. 
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yTs(P) = yh (9) +r {2t1-eni-cor-$) +2 [én ton +Z)-sin 2] 


, Nel 2c 
+> > aE oot sinn@-2n sin = cos ng 1} 


4 ast 4 (55) 


If the solutions for a particular power-law mean line are denoted 
by the superscript (a), such that the left hand sides of Eqs. 52 through 
55 become y . (x)/U,» Cp a) am, and ya) (9), then the solutions for 


a camber line represented by the polynomial of Eq. 37 are 


M 

y (Ve . yx, (56) 
M 

Cm= x Cm (2) (58) 
a=2 

sretoie 2 v8 9) (59) 
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SECTION III 


SOME NUMERICAL EVALUATIONS 


Several numerical evaluations are of interest. First is the 
determination of the impact of alternative quadrature rules and abscissae 
station distributions upon the accuracy of the integrals for a anat'°) 
given by Eqs. 31 and 34. Second is a determination of the relative 
contributions of the integrals for the forward, central, and aft regions 
of the airfoil as given by Eqs. 28 through 35. Finally, it is of interest 
to compare the aerodynamics as calculated by the Quadrature Method to that 
Calculated by the series-collocation approach used in the Power-Law 
Superposition Method. For the foregoing purposes, basic power-law camber 
lines normalized* to unity camber ratio, K= 1, are examined for exponents 
of a = 2, 4, and 8. Plots of some normalized camber lines are shown in 
Fig. 5. 


For the quadrature-rule and abscissae-station-distribution evaluations, 
comparisons are made for the trapezoidal rule and the quadratic rule of 
Appendix E in conjunction with dense (19 stations) and sparse (9 stations) 
distributions. Listings of the aforementioned abscissae distributions 


are presented in Table V**, 


There is no way of establishing the absolute accuracy of the previously 
mentioned alternatives for calculating uo) and Vc), However, "equivalent" 
collocation values of and Vc) are useful as reference quantities for 
assessing the "apparent" accuracy of the alternatives. A "collocation- 


equivalent value" is obtained by subtracting the quadrature-method 


* An aerodynamic property calculated for a normalized camber line may be 
applied to an arbitrary camber ratio for that camber line by simply 
multiplying the given property by the specified camber ratio. 


**Abscissae are listed only for the central region where numerical 
quadrature is employed. 
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leading- and trailing-edge contributions from the appropriate collocation 


solution quantity. The percentage relative difference in ae! is then 


defined as 
100(Hy!°)-(Hn!°)) squiy /(Hp!©)) equiy 


with a corresponding definition for gic): 


A tabulation of the percentage relative differences for ae and 
Mic) calculated by the alternative quadrature rules and abscissae station 
distributions for a = 2, 4, and 8 and e, values of .01, .10, and 4.00 is 
given in Table VI. It is seen in this table, that although there is some 
variation with "a" and Ness the most significant differences are due to 
the number of abscissae stations and the type of quadrature rule employed. 
It is apparent that the best results are obtained through use of the 
quadratic rule in conjunction with the dense station dist ibution, with 
the trapezoidal rule utilizing a dense distribution rating second best. 
The sparse distribution generally gives poorer results. Although the 
percentage relative differences for as with n # 0 are not shown in 
Table VI, those values presented for n = 9 are typical of those for n # 0. 
For a given computer capability and quadrature accuracy goal, Table VI 
provides a partial guide to the user in selecting the quadrature rule 
and station distribution most appropriate to his application. Another 
factor influencing the accuracy required in the numerical quadrature is 
the relative contribution of the quadrature calculation to the total 
value of 8. and Ac coefficients. In order to provide visibility on this, 
as well as visibility on the other components contributing to 8. and Ac» 
a computer printout showing the component breakdown is presented in Table 
VII*. Also shown, as a point of information, are values of Ac. calculated 
by both of the methods presented in the text. The Ac. value shown with 
the asterisk is calculated according to Eq. 17 and that without according 
to Eq. 25. The relative contributions of ne? to BF and go) to Ac. 
for the calculations shown in Table VII are summarized in Table VIII. 

It is observed in this table that, although there is some variation with 


won 


the exponent "a", the most significant variation is with c,, where it is 


F 
*The numeric quadrature calculations in Tables VII and IX employ the 
quadratic quadrature rule and the dense abscissae station distribution. 
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seen that the relative contribution of the central-region quadrature 
decreases considerably at the lower values of Cys Apparently, at low 


momentum coefficients, say c, = .10, the major contributions to 8. and 


J 
Ac. are from the trailing edge region; that is, from the terms containing 
the factor y'(1) in Eqs. 17 and 25, and the trailing-edge quadrature, 


Eqs. 32 and 35. 


Table VIII in conjunction with Table VI provides the user with 
additional guidance in the selection of a quadrature rule and abscissae 
station distribution. For os = .10 and a given accuracy requirement, 
the combined tables indicate a requirement less stringent than does 
Table VI alone. In some instances a relatively crude numerical quadrature 
calculation might be satisfactory. Consider, for example, that the 
interference lift coefficient (Acp = 4nB 5) for a camber line closely 
resembling an a = 2 power-law camber line is required to an accuracy of 
+2.5% or better at a ad 05. It is seen from Tables VI and VIII that a 


numerical quadrature using the trapezoidal rule and a sparse station 
distribution will satisfy this requirement (since -.3182 x 6.90 = 2.2%). 


As a final assessment of the quadrature method, collocation solutions 
for the same camber-line shapes and momentum coefficients considered in 
Table VII have been obtained. The resulting 8 and Ac. values are given 
in Table IX. The percent differences in Bo? Acp, and Ac. as calculated 
by the two methods are summarized in Table X. It is seen in this table 
that excellent agreement is achieved at the lower values of the momentum 
coefficient with the agreement becoming less favorable at higher momentum 
coefficients. There also is some degradation with increasingly higher 
Fourier harmonics. From examination of Table VI and VIII it can be seen 
that the momentum-coefficient degradation is likely due in part to the 
fact that the relative contribution of the numerical quadrature component 
to a total aerodynamic property increases more rapidly than does the 


quadrature accuracy for a fixed quadrature routine. 
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Although methodology is the main concern of this report, the 
interference lift- and pitching-moment coefficients for several power-law 
camber lines are presented in Fig. 6 as a matter of interest. As can be 
seen, for positive camber the interference of the jet sheet produces a 
lift decrement and moment increment, both of which increase with increasing 
jet-momentum coefficient and rearward movement of the maximum-camber 
location. The corresponding center of the interference lift is relatively 
insensitive to the camber-line shape and varies only slightly with the 
jet-momentum coefficient from a mid-chord position at low coefficient 
values to approximately the 46-percent-chord position at a jet-momentum 
coefficient of five. 
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SECTION IV 


CONCLUDING REMARKS 


A quadrature methed for calculating the incompressible-flow 
aerodynamics of an arbitrarily-cambered jet-flapped airfoil has been 
derived and evaluated. The principal advantage of the method is that 
it replaces the usual process of solving a set of simultaneous equations 
by numerical quadrature employing a table of influence functions. For a 
moderate-to-large number of simultaneous equations (say, nine or more) 
this results in a considerable simplification in the computation routine. 
As a consequence, relatively simple computational machinery, such as a 
pocket electronic calculator, can be used for computation instead of a 
digital computer. This facilitates convenient and fairly rapid analysis 
of prospective camber-line shapes. The method probably has its greatest 
utility for situations in which only a small or moderate number of cases 
are being examined. For extensive parametric studies, because of the 
sheer bulk of the computations, use of a digital computer probably is 


more desirable. 


Some numerical evaluations of the quadrature method were made by 
using aerodynamic properties calculated from a series-collocation method 


of solution as reference values. Factors examined were the effect of 


linear and quadratic approximating integrands and dense (19) and sparse 


(9) abscissae distributions on the apparent accuracy of the numerical 
quadrature, the relative contribution of the quadrature component to the 
total value of an aerodynamic coefficient, and the relative difference 
between aerodynamic properties calculated by the two different methods. 
The highest apparent accuracy in the numerical quadrature was obtained 
with the quadratic approximating integrand in conjunction with a dense 
station distribution, with the linear approximating integrand (trapezoidal 
rule) utilizing a dense distribution rating second best. The sparse 
distribution generally yielded poorer results. It was found that relative 
contribution of the quadrature component to a total aerodynamic coefficient 
increased with increasing values of the jet-momentum coefficient. For a 


given accuracy requirement on a total aerodynamic property this means that 


2k 
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the numerical quadrature accuracy at high speed (the regime of most likely 


application) can be less than that at low speed (higher c,). 


In general, for the cases examined, very good agreement was obtained 
between aerodynamic properties calculated by the quadrature and collocation 
methods (an average of .08% for the absolute value of the relative difference 
in the BF values) with the "apparent" error in the quadrature-method 
results increasing with increasing values of the jet-momentum coefficient. 
This partially verifies the necessity for increased numerical-quadrature 
accuracy as the jet-momentum coefficient increases. The tabulated results 
for the examples calculated in this report provide a guide to the user 
in selecting an appropriate abscissae distribution and numerical quadrature 


routine for a given application. 


Although methodology is the main concern of this report, a peripheral 
finding of interest is that for a positively cambered non-reflexed camber 
line the jet-sheet interference effect yields a lift decrement and moment 
increment both of which increase with increasing jet-momentum coefficient 
and rearward movement of the maximum camber location. The corresponding 
center of the interference lift is relatively insensitive to variations 
in camber-line shape and jet-momentum coefficient and is located approxi- 
mately at the mid-chord position. 
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APPENDIX A 


SOME RESULTS FROM UNBLOWN THIN-AIRFOIL THEORY 


It is the purpose of this appendix to present only some highlights 
of unblown thin-airfoil theory, including analytical expressions for the 


integrals required for polynomial camber lines. 


For a zero-thickness cambered airfoil at zero angle of attack the 
integral equation for the vorticity distribution 2) (x) is 


t +(o) 
e3 rine) 
y (x)= ae. f as dé 2) 


The solution (see, e.g., Ref. 4, 6, or 8) of Eq. A-~1 is 


Oh, 2 B(!8 PE é yl) 
Ue w\ x ) 0 (x- eves °° (A-2) 


for which 


of) =-27a'g) 


cl = oy +a, (A-k) 
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In applying Eqs. A-2, A-5, and A-6 to arbitrary camber lines, 
numerical quadrature is usually employed. For bounded values of y'(0) 
and y'(1), evaluation of the foregoing integruls gives no difficulty at 
the leading- and trailing-edge singularities, although approximate 
analytical treatment paralleling that of Appendix D likely will be 
necessary for small regions near these extremities. For Eqs. A-5 and 
A-6, this subject is discussed in more detail in Ref. 7. Special treat- 
ment in the region of € near x also will be necessitated. Studies of the 
numerical evaluation of Eq. A-2, or alternative forms thereof, may be found 
in Refs. 5, 7, and 9. 


The principal purpose of this appendix is to present some pertinent 
relations for polynomial camber lines. For polynomial camber lines of 


the form of Eq. 36, the following basic integral repeatedly occurs in 
Eqs. A-2, A-5, and A-6 


I 
ed ee 
ce Sresy areaths 


M2 O12; anna, GEXSI (A-7) 
where the integral is the Cauchy principal value. 


The integral is evaluated by 


Ten (X) = XLpm (x) + Ln 44 (0) (a-8) 
where 
Ig(x) = 0 (A-9) 
T, (0) =-7 
Im+) (0) =-1-3-5----_ (2m-1)w/2™m! (m#0) (A-10) 


The integral of Eq. A-10 is taken from integral 212.4 in Ref. 10. The 
integrals given by Eqs. A-8 through A-10 are analogous (but not identical) 
to integrals 6 through 13 in the Appendix B of Ref. 8. 
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APPENDIX B 


THE FUNDAMENTAL JET-FLAPPED AIRFOIL SOLUTIONS 


The jet-flapped airfoil flows considered to be "fundamental" for the 
purposes of this report are the flows for the regularly blown flat plate 
(Ref. 4), the singularly blown flat plate (Ref. 4), and the regularly blow, 
mechanically flapped airfoil* (Ref. 6) (see Fig. 1). Since most of the 
mathematical relations involved in the latter two cases are employed or 
implied in the methodology of this report, it is of interest to summarize 
them in this appendix. 


SINGULARLY BLOWN FLAT PLATE 


This case is illustrated in Fig. lb. The pertinent relations as 
given or implied by Ref. 4** are as follows 


‘ N-I 
Yq (0) X)= 2g, x2? [ fn (1x) + 2AQX/U1 +X) + X AnX" | (B-1) 
n= 
Cer(cj) = 477 Ay (B-2) 
N-I 
Cmr(cj)= ~ (Cj +14) - X Antn (B-3) 
n= 


* Called a "jet-augmented flap" in Ref. 6. 


**The relation for the chordwise vorticity distribution in Ref. 4, Eq. 107, 
contains typographical omissions. The second and third terms on the right- 
hand side of Eq. 107 should be multiplied by x73/2 and the term 
2af1-x)/x]1/2 should be added to the equation. 
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| 
lg =-(4rif x! fn (i-x) dx=(16/m)(t-n 2) 


| 
lo = af x'/2 (1-./T=x ) dx= a(i- 5) 


Vi 
In= aff «2X ax (n=1,2,3,---) 


where A) = A,(cs)s X is defined in the list of symbols, and 
1, = 4(4-2-K)) 


In = a[-(2n+1y/(2n-1) + Ky. | 
with a given by the recursion formulae 

K,=4-7 

Ky = /(2n-th-Ky-| 


The AL coefficients are given by 


N-I 
L_ (Oma + ABmn) An =m + Adm 
n=0 
where 
Amo = SING m (n=O) 
Omn = (1 + Cos@m) SINGm (n>1) 


9 
Din = 4(cos ng+2n tan — sin ngm)/(4n®=1) 


(B-4) 


(B-5) 


(B-6) 


(B-7) 


(B-8) 


(B-9) 


(B-10) 
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dey = (8/7) {| 806(%m/2) | 4n[ tant 9 /4)]-Anl tant aim /2)] } (Bel1) 


Qm=mmr/N m=0,1,2,----,N-! (B-12) 


The equation for the trailing streamline is given by Eq. 123 in 
Ref. 6. 


REGULARLY BLOWN, MECHANICALLY FLAPPED AIRFOIL 


This case is illustrated in Fig. le. The pertinent relations as 


given or implied in Ref. 6 are as follows: 


% (cj,x,€)= % SNe) +», (4) (Gx, €) (B-13) 


(0) 2u sin | 

y (x,€)2 =" dy ton ZS +h . 
8 ‘~ sin 2X (B-14) 
7 
N-I 
% A (ej,x,€) = 20, £8 [200K 14X) +E onX"| (B+15) 
Dnt Dale), €) ; 
xe (i/2itcosO); = (1/21 +cosyx) (B-16) 
Cas (¢j,€) = cep! (€) +47Do (B-17) : 

(0) No 

emg (Cj, €) = omg (€)-€ej-E Dn In (B-18) | 
| 
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where I is given by Eq. B-6. The Do coefficients are given by 


N-l 
2 (ama * ABmnn) Dn (cj, €)= tm(€) (B-19) 


where m = 0, 1, 2, ----, N-1, . and a are given by Eqs. B-10, and 


2 Pm 4 bd “1 tan 
fm(€) = —& tan — ~sec —@ tan? |——&= 
x Pon in ae sina (B-20) 


As discussed in the main body of the report, correct implementation 
of the quadrature method requires that the aerodynamic derivatives Cp (c, 58), 
6 


ch (c, 6) and Vgc, x58) for the blown mechanically flapped airfoil reduce 
6 


to the corresponding derivatives with respect to a for the regularly blown 
flat plate and with respect to Tt for the singularly blown flat plate for 
— = 0 and € = 1 respectively. 


That is, the following must be satisfied: 


Chg (c},0) = Cgy (¢}) 

Sms (cj 10) = Cma(cj) 
;. = : -21 

% (c),*,0) % (ej ,x) (B-21) 

Chg (Cj,!) = cg (cj) 

Sma (cj, = Cm _(C}) 


% (cj xe y (¢},x) (B-22) 
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For — = 0, Eq. B-16 yields x = 7, and f (0) of Eqs. B-20 becomes 
e, Where e. (defined in Ref. 4) is the correct right-hand side for Eas. 
B-19 to yield the numerical results D,(¢, 20) = Bale,)s where the B.'s 
are the coefficient of the regularly blown flat solution (Ref. 4). Hence 
in this case the proper limit is obtained. 


We now consider the limit — = 1 (X = 0). The proper limit for f,(2) 
is obtained by noting that f,(E) = -2, for which f,(1) =-2. The result 
f,(1) = 0, obtained by first substituting — = 1 and then % = O in Eqs. 
B-20, is incorrect. The foregoing choice for f,(2) also may be verified 
by checking the trailing-edge boundary condition. From Eqs. 17 and 23 
in Ref. 7 it may be shown that this boundary condition is 


N=! 
D_ Dp (ej 1971 4n®-1)=-1/20 
n=O 


Substitution of f,(1) = -2 in the first of Eqs. B-19 satisfies the boundary 
condition, whereas f,(1) = 0 does not. We also have f (2) = 0 for m #0. 
The coefficients D,(e, +2) are then determined by 


FY Maiti. } 
dig [Gen * AP mn) OnE} Y= 1 o me1,2,3,---- 


Also for — = 1, ne (x,1) = 0; then 


v (a) 


% 


(¢),x,1)2 2U, 0°" {af ngt¢},1)] XAI +X0+5 [D4(ej0]X"} (B-2h) 


where the inverted circumflex indicates that the vorticity distribution 

is obtained by taking — = 1 in the regularly blown, mechanically flapped 
airfoil solution of Spence (Ref. 6). From a visual comparison of Eqs. B-9 
and B-23 and Eqs. B-l and B-2h, it is not possible to determine whether 
the conditions of Eqs. B-22 are satisfied. However, extensive numerical 
calculations for N = 3 and 9 reveal that they are not. Because of this 
discrepancy, quantities obtained by taking —€ = 1 in the Spence theory will 
be denoted by an inverted circumflex to distinguish them from the correct 
result; for example, Cp(6), 1), Emp (C4 o!), and Yalejnnt 
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THE INCREMENTAL COEFFICIENTS AD AND AS 


The incremental coefficients AD, and AS are defined as 


ADq(¢j) ® Aq (cj) -Dp(¢js!) (B-25) 


(B-26) 


N-I 
AS ® Sy~Sp=Z =ADnIn 
n=0 


where Sys S)» and I, are defined by Eqs. 23, 24, and B-5 and B-6. These 
incremental coefficients can be determined by separately evaluating A, 


and Do from Eqs. B-9 and B-19, or, in view of the linearity of the problem, 
by solving 


N-I 
2+c7+Ad m=0 
{damn + NOmn) BOn = { cm +\dm i (B-27) 


for AD. 
n 
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APPENDIX C 


DETERMINATION OF Dt AND SS 
The second derivative of 2. with respect to &, pe is determined 


from the set of simultaneous equations resulting from taking the second 
derivative with respect to € of both sides of Eq. B-19 yielding 


"5, (mn XP Dn = im| €01-€) |? 


(c-1) 
where 
| @ 213-26) 2(1-€) 
len * ht [ue ~« 
| me tone | 2) eee) * (i -EeoeP DE me 


ha Kc oS et mA INRA Rar ee enn ae 


} The right-hand side of Eqs. C-l has singularities at §€ = 0 andl. In 
order to avoid numerical difficulties associated with these singularities 
| it is convenient to introduce the coefficients En defined by 


Dy we, | €t1-€)] > (c-3) 


Equations C-1 then yield the set of simultaneous equations 


he ca ek UREN 


N-I 
Ro (mnt Amn) En = he (c-4) 


which are solved for the E.'s: 


The relations for so: given t_- Eq. 27, may be written as 


Sp: s[Eu-e] 


N-| 


= E,!I 
s 2X nln (c-6) 
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APPENDIX D 


TREATMENT OF THE LEADING- AND TRAILING-EDGE 
SINGULARITIES FOR THE QUADRATURE METHODS 


LEADING EDGE 


Evaluation of the following integrals is required: 


tines [ve Dn(€) d€ 


Sark E88 (Erde 


where a << 1 and 
Dn (€)= En(€)[€(1-€)] 
spi€) = s(€)[€c1-€)] 


with EC) and s(&) given by Eqs. C-4 and C-6. 


Numerical calcwlation shows that E(e) and s(&) are well behaved in 
the vicinity of € = 0 and therefore can be approximately represented by a 
quadratic equation in — in that region. It is assumed also that y(&) is 
approximately a quadratic equation in the same region. Since y(0) = E (0) = 
s(0) = 0, we may write 


En(€)= a, € +a,€* (D-5) 
s(€)=b,€ +b,€" (D-6) 


yl€)=c,€ +e2€* (D-7) 
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The constants B)> Bp» dy» and db, are determined by passing Eqs. D-5 and 
D-6 through the known values of E (Eq) » E(E)s s(E). and s(E,), where 
Oo< Eo < E- The constants cy and c, are determined by constraining Eq. 
D-7 to yield the known values of y'(0) and y(E,). Determining the 
constants as specified, substituting Eqs. D-5 and D-7 in Eq. D-l, and 
substituting Eqs. D-6 and D-7 in Eq. D~2 yields 


(HaLeE ={ LoEntEd +h, En(€)| WE) +{ MoE n(Eo) +m, En(€,)] y'(0) (D-8) 


Mie + [Lo sé) +4, 9(€i)] v(€,)+[ mos( Eq) +m, 8(€)]y'(0) (D-9) 
where 
£. = 2381-4 - 7 Kotys 

. oS 1E,-Eo) ” €7(€,-&o) (D-10) 

mae Jeb -2u3Ei tq no 2 +ylEg+€,)-d 
€0€s(€,- Eo) "tee (D-12) 

i = . 

ine’ €"[€01-€] 7 a€ (D-12) 


J = 2/8, 70-8) ( for n=l) 


Jn= Sn-} = In-2 ( for n>1) 


where 


jo =-sin' (1-2€,) +972 
in ®- JE-&) +1972 
ie? -4[ 6 +(3/2)| VE-E) + 31978 (D-13) 


; 
. 
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In the present application bo and a were taken as 55 = 0.025 and 
B = 0.050. These values yield 


fo = 2.1337 £,= 0.82274 


Mo = 0.13977 m, = -0.0082727 (D-14) 


Numerical values of BA(c5 58) > BA(c,.8)), s(c5sE). and s(c, 8) for the 
above values of bo and a are given in Table III. 


TRAILING EDGE 


Evaluation of the following integrals is required: 


(Hn) TE -f y(€)0n (€)d€ 
€2 


res f. viersy (eae 
€e 


where (1 - Eo) << 1 and DY and s5 are defined by Eqs. D-3 and D-4 


respectively. 


Analogous to the leading-edge development, E,(6)> s(&), and y(&) are 
approximately represented by quadratic equations in the vicinity of the 
trailing edge. 


Since y(1) = 0, we may write 


En(€) = En (t) + a, (1-€) + agtt-€)? (D-17) 
s(€) = (1) +b, (1-€) + ba (t-€)? (D-18) 


yWie)= c,(1-E) +eg(1-€)? (D-19) 


The constants a> 5s bi» and b, are determined by passing Eqs. D-17 and 
D-18 through the known values of BE (E)> EA(E3)> s(E,), and s(&,), where 


Eo < b3 <1. The constants c, and c, are determined by constraining 
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Eq. D-19 to yield the known values y(E,) and y'(1). Determining the 
constants as specified, substituting Eqs. D-17 and D-18 in Eq. D-15, and 
Eqs. D-17 and D-19 in D-16 yields 


(Hn) TE = [4 Enl Ep) + 4s En (Es) + 24 En( 0] y(&2) 


+ [ moEn( &2) +m3Eq( Es) + mgEn tt] y'(1) (D-20) 


Are s [ 42812) +h; s(€3)+4,3(0)| y(&>) 


+ [mzs(€_)+mgs(Es) +mas(l)] (1) (D-21) 
where 
a ~U3(1-§s) 404 
2-82) 5(E5- Ee) 
$3 Jslt-&a)-uUg 
5 1-&)? = Es Ea- Ep) 


. dg 1-2) 1-&s)-U3(2-€ 2-3) tug 
(1-€)°U1-§3) (D-22) 


7 U2 (I~ Ea1-&3)-J3(2-E2-Es) +4 
(1-€2)"(€3-€2) 


_ Yall &e)*+2 3 (1-E 2) ~via 


: i (1G git- Sg y-G) 


. = Jj (1 - e)7(1-€ 3) +02 (1-Ep 3-€ 2-265) -d 33-28 0-E5) +g 


ia (1-€)°(1-€5) (D-23) 


(D-2h) 


Un® i 1" €1-E)] 7 aE 
2 


H 
f 
i 
4 
# 
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Ui = 2A (1-&,)7€, 
U2 =u) -jo 
Js 20, -2i0 ti 


42d - Slot 3ii ie 


jo= 7/2 + sin'(1-2€>) 
i = SEglt-Se) + 10/2 
ia = (4 E243) VEq(I-&) + (3/8)i0 


In the present application Eo and B3 were taken as Eo = 0.975 and 
E, = 0.9875. These values yield 


L, = 0.82274; A, = 2.1337; , = 0.16545 


ma = 0.0082727; m= -0.13977; mg=-0.17124 (D-26) 


Numerical values of Eile, 56) Bic, 563), Bi(e,51), s(c, 65), s(c, +63)» 


and s(c, 51) for the above values of E> and b3 are given in Table III. 


EEE SII mY 
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APPENDIX E 


NUMERICAL QUADRATURE ALGORITHM 


A numerical quadrature formula for an arbitrary number of unequally 
spaced abscissae is developed in this appendix. 


Referring to Fig. 4, the integral under consideration is 


ansxo)® f™ (xx 


(E-1) 
This is approximated by 
N-| 
L(xy X0)= DL Alp 4 (E-2) 
n=0 
where 
Xn+l 
Aln+) = f(x) dx 
walt. rs (E-3) 
and 
f(x) = fq + D(x-xq) + 6 x-xq)? (E-h) 
yielding 
Alyn 4] = fl Xn 4i-%p) + (D/2)(Xp4{-Xq)? + (6/3) Xn +1 -%,)° (E-5) 


Two of the collocation points for the approximating parabola of Eq. E-4 
are taken to be the ordinates of the sides of the elemental area AT ae 
The form of Eq. E-4 assures collocation at x? whereas collocation at 

X41 will be achieved by proper determination of the coefficients b and c. 
A third collocation point is needed, say at Xi» to ensure the determinacy 


of b and c. These coefficients are determined such that 


B= (XH 44, %j) , C= (Ky 41, Xj) (E-6) 
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where 


‘ { n+2 “Forward collocation” } 
is 


n-l "Backward collocation” (E-7) 


With this scheme, there are two options available in applying the composite 
quadrature formula of Eq. E-2. For the first option backward collocation is 
applied to all the elements except the first, at which forward collocation 

is used. For the second option forward collocation is applied to all the 
elements except the last, at which backward collocation is used. The first 
option is employed in the present application, since in an automatic hand 
computation errors of omission are less apt to occur if exceptional operations 
are performed at the beginning of the computation rather than at the end. 


Collocation at x» x and x; yields the result 


n+1? 


N-l = — 
I(xy.xo)= 5 nar Lory | ae Picea (frain fa) 


n=0 = 8 (xj-Xp+1) 
sj] 
cx fr eee dee BER S 
| caoltgecte (fi- Fn) (E~8) 
with 
i=2 for n=O 


izn-l for n>O 


A closely related quadrature formula is given by Davis and Rabinowitz 
on page 48 of Ref. 11. The Davis-Rabinowitz formula uses forward collocation 
for the first element, backward collocation for the last element, and an 
average of forward and backward collocation (called "overlapping parabolas" 
in Ref. 11) for the central elements. The method of Eq. E-8 entails 
considerably less computational labor than the method of overlapping 
parabolas, and is to be preferred from that standpoint. Also, it is 
demonstrated in the numerical evaluations section of this report that the 
present method is more than adequate for most applications. Additionally, 
it was surprising to discover that for the examples on page 50 of Ref. 11 
the present method yields errors: smaller than the overlapping-parabolas 
and the cubic-~interpolation methods in that reference. Depending upon 
the case, errors of 10 to 50 percent smaller were obtained by the present 
method. 
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It is not the intention here to explore the merits of various 
numerical quadrature methods in depth. It is believed that the present 
method represents a reasonable compromise between accuracy and complexity 
suitable for the purpose of this report. The user may, of course, employ 


any other method of his choice. 


The selection of the quadrature formula of this appendix was governed 
by considerations of use in automatic hand computations. As a matter of 
convenience and to ensure a consistent basis of comparison for hand and 
machine computations, however, the same formula is also used in the 


digital computer program. 


a 
: 
i 
: 


ee — 
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APPENDIX F 


NUMERICAL QUADRATURE PROGRAMS FOR THE HP-25 AND SR-56 CALCULATORS 


The algorithm presented in Appendix E is implemented in this appendix 
with program listings for the Hewlett-Packard HP-25 and Texas Instruments 
SR-56 calculators. These calculators typify the commonly used Reverse- 
Polish and algebraic logic systems. 


In the run instructions given below for the aforementioned programs, 
the following symbology is used. A nonunderlined term is the symbolic 
representation of numeric data to be input. An underlined term is an 
exact or abbreviated representation of the symbol on the calculator key 


to be pressed. 


HEWLETT-PACKARD HP=25 


Program Listing 

1 x 11 ST0-0 21 RCL 1 3) Ws 

2 sto 6 12 RCL1 22 RCL 2 32 x 42 x 

3 RCL 5 13. - 23 «+ 33 + 43 ST0+7 
4 - 14 x2y 24 + 34 RCL 5 4A RCL 2 
5 xey 15 + 25 RCL1 35 STO 4 45 $17 1 
6 sw O 16 2 26 RCL 3 36 6 46 RCL6 | 
7 ST 3 17 + 27 - 37 = 4? ST 5 
8 RCL1 18 x 28 + 38 + 48 RCL3 
9 - 19 RCLO 29 RCL 5 39 RCLO 49 ST 2 
10 RCL 2 20 g(x“) 30 RCL 4 6 


Run_ Instructions 

To initialize, key: 

1. X2_ ST 1, ¥os 4 Bs X» ST 4 

2. Xo» STO 2, Yor ts So» X» ST 5 

To run, key: 

3. ™%» aa) yy» 4, G1» R/S (xy is displayed) 

4. Repeat run step 3 for XQ» Xv - = = =» Xy and corresponding values of y and 


The current abscissa input, x,, is displayed each time run step 3 is 
performed. This serves as an orientation aid in inputing the data. 


5. Press RCL 7, read I(xy, x). 
yy 
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TEXAS INSTRUMEN'S SR-56 


Program Listing 


R/S 
STO 
3 

SUM 


e 
e 
sad 


cr 


WOMBIYMMAWNEWNH—O 


a 


Bt 4oRgnk-wg! Fz 
e 
o 


oe 


28 
29 


Run Instructions 


At ORZtrngt+ von ing 
i ad 

A AtORnN ing at ~og 
e 


BSSYRREYVBYS SES ES 


ES 
7x 


To initialize, key: 


1. X25 ST 2, ST 8, yo» Xs 2» =» ST 4 


Xo» STO 2, +/-» SUM 2, STO 9s Yo» X» Bo» =» ST 5 

run, key: 

X1» R/S, yy» R/S» ey» R/S (x, is displayed) 

Repeat run step 3 for X20 Xr - eo me XN and corresponding values of y ard g. 


The current abscissa input, x,, is displayed each time run step 3 is 
perfcrmed. This serves as an orientation aid in inputing the data. 


Press RCL 0, read I(xy,Xq). 
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APPENDIX G 
COMPUTER PROGRAM SYMBOLS * 


Sinee the symbols used in the computer printout and in the program 
internal logic are not necessarily related or compatible, separate PRINTOUT 


and FORTRAN symbol listings are given in this appendix. 


PRINTOUT SYMBOLS USED IN TABLES I, II, III, IV, AND IX 


Symbo 1 Definition 
A power-law camber line exponent a, see analysis symbols 
A(0),A(1), ete. see analysis A, 
BETA(O),BETA(1),etc. see analysis B,, 
Cc denotes the central portion of the airfoil 
05 <x < .975 over which a numerical integration 
is performed (Table VII), see analysis ( fe) 
CJ see analysis c 5 
C(0),C(1),ete. see analysis C,, 
DCM see analysis Ac, 
DD(0) ,DD(1) ,ete. see analysis AD, 
DDP(0),DDP(1),etec. read as "dee double prime," see D. 


DELTA D(N) see analysis SD, 


“Generally, the computer symbols have been defined in terms of the analysis 
symbols for which definitions are given at the beginning of the report. 


' 
{ 
t 
; 
5 
k 
4 
: 
§ 
' 
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PRINTCc? SYMBOLS (Contd) 


DELTA S 
D(N) 

Ds 
E(0),E(1),ete. 
ET(M) 

FT(M) 
H(0),H(1),ete. 
H(N) 

HC 

I(N) 

ISCRPT 

KAPPA 

Ics 

LE 


RA 
SDP 


SUMI(N)A(N) 


SUMI(N)C(N) 


Te 


TOT 
XI 


_ anaes 


see analysis AS 
see analysis D,, 
see analysis AS 
see analysis E, 
see analysis e,, 
see analysis f = 
see analysis H, 
see analysis H, 
see analysis ale? 
see analysis Lh 

see analysis c 

see analysis K 

read as “lower case sum," see analysis s 
denotes the airfoil leading-edge region, 
O0<x< .05 (Table VII) 

see analysis r, 
read as "S double prime," see analysis S" 
N-1 

2% TyAn? see analysis symbols 


N-1 

2 InCn» see analysis symbols 

n= 

denotes the airfoil trailing-edge region, 
975 < 0 < 10 (Table VII) 

denetes a total value 


see analysis & 
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PRINTOUT SYMBOLS (Contd) 


YPDD(0),YPDD(1),etc. read as "y prime delta dee zero," etc. 
i.e., y'(1)ADp, see Eq. 11 
YPILDS equals y'(1)(Ip+AS), see Eq. 25 


FORTRAN SYMBOLS USED IN PROGRAM LOGIC AND PRINTOUT 


Symbol Definition 
A exponent "a" for basis power-law camber line 
A(M,N) see analysis a 
B(M,N) see analysis b, 
BETA see analysis B,, 
BETAC program logic flag for collocation solution BETA 
calculations 
BETAQ program logic flag for quadrature method BETA 
calculations 
N-1 
BISUM xX BT,’ See analysis symbols 
Bl 4/ (un? - 1), see analysis symbols 
same as C in PRINTOUT SYMBOLS 
a ~*~ Ab 
mn mn 
cc same as C 
CJ see analysis c j 
CJVSE see program input definitions 
CMPHAI see program input definitions 
CONST array of constants used in H and ISCRPT calculations 
DCM see analysis Ac,, 
DDP E( )/DIVSR( ) 
DDT see EA below 
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FORTRAN SYMBOLS (Contd) 


DIVSR divisor value array 
divisor value input 
denotes ©, + df, (in analysis symbols) before simu1- 
taneous equations are solved and the Fourier coefficients 
A,» Dy» Bp» etc. after the simultaneous equations are 
Solved 

EA( LSE) alternative symbol for Fourier coefficients, A,, Dy» 
Ey» ete. 

EB( LSE) alternative symbol for e, + Af, 

EPS (LE, ISE) see analysis e,, 


EPSI( LSE) see analysis en 
ET(M) see analysis en 


FIRST program logic flag identifying the first pass through 
the program 


program logic flag related to BETAQ 

n, n=1, 2, ---, N-1 in analysis symbols 
see analysis fn 

see analysis tn 


see analysis HL 


see analysis H (ce) 
incremental HC value 
array containing all of the output page headings 


see program input definitions 

see program input definitions 

see program input definitions 

see analysis (Ay) rp 

see analysis (Hy) ass 

error flag returned from simultaneous equation solver 
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FORTRAN SYMBOLS (Contd) 


IINTRR incremental ISPTC value 

IL see analysis ly 

IN see analysis a 

ISCRPT see analysis 

ISP see analysis g'° 

ISPTLE see analysis So 

ISPTTE see analysis Son 

KAPPA see analysis kK 

KHLE bracketed factor of y'(0) in Eq. D-8 

KHTE bracketed factor of y'(1) in Eq. D-20 
KTIE bracketed factor of y(&;) in Eq. D-8 

KTTIE bracketed factor of y(€,) in Eq. D-20 
LAMBDA see analysis 

LE a program loop variable related to XI or A 
LL a program loop variable related to LAMBDA 
LSE a program loop variable related to the number of 


simultaneous equations 


ISEA see ISE above 

NE the number of XI's or A's 

NEM3 NE-3 

NEM4 NE-4 

NL the number of LAMBDA values 

NSE the number of simultaneous equations, see analysis N 
OPTION see program input definitions 

PHI see analysis 9, 
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FORTRAN SYMBOLS (Contd) 


YPILDS 


nm (to 13 decimal places) 

power-law camber line exponent a, see analysis a 
see analysis r, 
same as in PRINTOUT SYMBOLS 
bracketed factor of y'(0) in Eq. D-9 
bracketed factor of y"(1) in Eq. D-21 
bracketed factor of y(€;) in Eq. D-9 
bracketed factor of y(€,) in Eq. D-21 
sinn@,, see analysis symbols 


N-1 
2% TE, where E, may be any one of the Fourier 


coefficients, Ay, By» Cy» Dy» etc. 

SUMIE( )/DIVSR( ) 

same as in PRINTOUT SYMBOLS 

same as in PRINTOUT SYMBOLS 

XI or A value array 

XI or A input value 

see analysis € 

array of airfoil camber-line displacement coordinates 
y'(1)4D, 

slope of the camber Line at the airfoil leading edge 
slope of the camber line at the airfoil trailing edge 


equals y'(1)(Iy + 4S), see Eq. 25 
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APPENDTX H 
DIGITAL-COMPUTER PROGRAM 


A listing of the program written to calculate the parameters discussed 
in this report is given in this appendix. The program is written in. 
FORTRAN IV language and consists of a main program and two subroutines. 

The program MAIN reads the input and performs the calculations. The 
subroutine WRITE handles all of the output printing. The subroutine SIEQ 
solves a set of simultaneous equations*. Since most computer systems have a 
simultaneous-equations library routine, a listing of SIEQ is not given. A 


discussion of the output options and input parameters follows. 
PROGRAM OUTPUT OPTIONS 


(1) Expanded output (A(M,N) and B(M,N) along with some intermediate calcula- 
tions are printed). This option is principally for diagnostic purposes. 
There is no example in the report. The option is implemented by setting 
OPTION = .TRUE.. All other logic flags must equal the default values. 

(2) Standard output (see Tables I, II, and IV). This is the default output. 
All logic flags must equal the default values. 

(3) Output similar to standard output but with XI varying at a constant value 
of CJ (see Table III). This option is implemented by setting CJVSE = 
eFAISE.. All other logic flags must equal the default values. 


(4) Output as in option (3) with values of the airfoil ordinates and 


*The reader should note that in the argument list for subroutine SIEQ the 
array, E, containing the column matrix is redefined during the subroutine 


execution and is used to return the Fourier coefficients. 


eG ae eit ee i eee ee in 
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the parameters H and ISCRPT computed and printed (no example in report). 
This option is implemented by setting CJVSE = .FALSE. and CMPHAI = .TRUE.. 
All other logic flags must equal the default values. 

(5) Airfoil ordinates and the parameters H, ISCRPT, YPDD, and BETA for the 
quadrature method computed and printed (Table VII is an example of this 
type of printout except that the airfoil coordinates have been omitted). 
This option is implemented by setting CMPHAI = .TRUE. and BETAQ = .TRUE.. 
All other logic flags must equal the default values. For this output 
option the first set of HEADI and $PARAM3 cards must be for the case of 
4D,» where AD,, is defined in Eq. B-27. 

(6) BETA output for the collocation solution in the same format as option (2). 
This option is implemented by setting BETAC = .TRUE.. All other logic 


flags must equal the default values. 
PROGRAM INPUT 


For ease in use, the NAMELIST input format is employed to input all 
numerical values. The input namelists and heading cards must be in the 
following order: 

The NAMELIST Group-Name $PARAM1 

NSE number of simultaneous equations 

NL number of lambda values (maximum of 30) 

LAMBDA (1) list of lambda values 

Various output options are implemented by the following logic flags when used 
in accordance with the PROGRAM OUTPUT OPTIONS described above. 


CJVSE «TRUE. = output in the form of one XI for range of CJ's (default) 


-FAISE. - output in the form of one CJ for range of XI's 
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CMPHAI eTRUE. - H's and ISCRPT's are calculated, requires $PARAM2 


-FALSE. - H's and ISCRPT's are not calculated, $PARAM2 must 
not be present (default) 
BETAC eTRUE. - Beta values for collocation solution are calculated 
-FALSE. - Beta values for collocation solution are not 
calculated (default) 
BETAQ eTRUE. - Beta values for quadrature method are calculated, 
requires CMPHAI = .TRUE. 
-FALSE. - Beta values for quadrature method are not calculated 
(default) 
OPTION TRUE. - A(M,N) and B(M,N) arrays are printed along with some 
intermediate calculations 
-FAISE. - A(M,N) and B(M,N) arrays and intermediate calculations 
are not printed (default) 
The Heading Card HEADH 
This card contains the heading used for the printout of the airfoil 
coordinates and the parameters H and ISCRPT (maximum of 80 characters centered 
on card) 
The NAMELIST Group-Name $PARAM2 
Y(1) airfoil camber-line ordinates used with XI's, beginning with 
XI = .05 and ending with XI = .975 (maximum of 30 values) 


YPLE slope of the camber line at the airfoil leading edge 
YPTE slope of the camber line at the airfoil trailing edge 
PICLA power-law camber-line value 
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The Heading Card HEADE 
This card contains two headings; one for the Fourier-coefficient array 


(maximum of 5 characters right justified in columns 1-5) and a second one 
for the summation array (maximum of 15 characters in colums 11-25 centered 
about column 20) 

The Heading Cards HEADI 

Two mandatory cards are required to provide the heading for the output pages 
(maximum of 136 characters; 80 columns of the first card and 56 columns of 
the second card). These cards must precede each occurrence of the $PARAM3 


namelist. 


The NAMELIST Group-Name_ $PARAM3 


VARI XI or the parameter A for the power-law camber line 
DIVSRI divisor value for corresponding VARI value 

EPSI(1) represents the analysis symbol e,, 

FTI(1) represents the analysis symbol f,, 


The HEADI and $PARAM3 namelist cards are repeated for each XI or A. 


PROGRAM LISTING 


The program listing is given on the following pages. 
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PROGRAM MAIN 


1 PROGRAM MAINCINPUT ,OUTPUT,TAPES=INPUT yTAPEG=OUTPUT, PUNCH) 
REAL LAMBDA,KAPPA,IL,IN 
REAL ISPTC,IINTR, ISPTLE,ISPTTE,ISCRPT,KTLEsKHLE, KITE, KHTE 
COMMONZALL/ A959) BO 9,9) C1949) 5001959) ,E (9) 5EB19) CU (9) 
5 19EA03059) CPS (3059), FT (3059) »HEAD(30 517) y DOP (305 9) » SUMIE( 30) ,SDP(3 
20) s0IVSR(30), JAR(30) sCUdSE,LAMBDA(30) yNSE gNLyNEs OPTION,FIRST 
3,HC (9) sISPTC,Y(30) sCMPHATI sHLE (9) yHTE (9) ,H(9), ISPTLE, ISPTTE,ISCRPT 
4yHEADE (5) »yPLCLA,HEADH (10) »BETAQ,BETAC,BETA(9) y YPOD(9) y FLAG, 0CM 
5, KAPPA,RA,0S(30) ,OCMASK, YPILDS 
10 INTEGER PLCLA 
OIMENSION EPSICY9) »FTI(9) y INC9) »HEADI (17) y CONST (10) OO T(30,9) 
LOGICAL OPTION,CUVSE,FIRST,CMPHAI » BE TAQ, BETAC,FLAG 
NAMELIST /PARAML/NSE yNLy OPTION sL AMBDA,CUISE sCMPHATs BETAQs BETAC 
NAMELIST /PARAM2/Y,YPLE,YPTE,PLCLA 
15 NAMELIST /PARAM3/DIVSRI,EPSI,FTI, VARI,KAPFA 
DATA CONST/2.1337 50082274 50413977 »-060082727,1012293,5 
R2 0952695 042377950. 00594475,-0.0975962, -0012077/ 
DATA IN/1671081551613274150. 401184, 0619822 450.11 6654, 0.076405, 
A0 053769, 0. 039834,0.030666/ 
20 c 
Crt tte INITIALIZE LOGIC FLAGS 
Cc 
FLAG=eFALSE. 
OPTION=eFALS=. 
25 CJUVSe=.TRUc. 
CMPHAI=.FALSL. 
BETAQ=.FALSE. 
FIRST=eTRUZe 


30 KAPPA=1.£60 
Cc 
crs. OEFINITIONS 
Cc . 
NE=0 | 
35 PI=3.1415926535896 


IL=(16-/PI) *(1e-ALOG(2.)) 


40 Uteeee RcAD INPUT DATA 


READ (SsPARAM1) 
IF(BETAQ) FLAG=.TRUE, 
IF(CMPHAT) READ(5,1100) (HEAOH(I) »1=1,10) 
45 IF(CMPHAI) READ(5,PARAM2) 
READ(5,1110) (HEADE(I),1=1,5) 
10 READ(5,1100) (HEADICI) ,1=1,17) 
IF(EOF(5)) 40,12 


12 NE=NE+1 
50 READ (5yPARAM3) 
c 7 
Grr rte OcFINE INPUT BATA ARRAYS | 
c ; 


VAR(NED=VARI 
55 DIVSRENE) =OIVSRI 
OO 20 LSE=1,yNSc 
EPS(NE ,LSE) =EPST(LSE) 


BEST AVAILABLE COP 


56 
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PROGRAM MAIN 


60 


65 


70 


75 


60 


85 


90 


95 


105 


110 


20 FT(NE,LSE) =FTIC(LSE) 
DO 30 K=1,17 

30 HEAD (Ney K) =HEADI(K) 
IF(FLAG) GO TO 40 

35 IF ttor(5)) 40,10 

40 CONTINUE 


c 
Clever CALCULATE A(MyN)*S AND B(M,N)*°S 
c 
00 90 M=1,NSE 
PHI = ((M-1)*PI)/NSE 
A(M,1) = SIN(PHI) 
B(Ms1) = —4e 
00 90 N=2,NSE 
FN = Nel 
SN = SINCFN*PHI) 
A(MyN) = (ie + COS(PHI)) *SN 
Bi = 4e/(4e*FN*FN - 10) 
SO B(MyN) = Bi*(COSCFN*PHI)D + 26 *FN*(TANCPHI/20)*SN)) 
c 
CESete TEST FOR EXTENDED QUTPUT OPTION 
IF(OPTION) GO TO 95 
c 
Gee see TEST FOR OUTPUT IN FORM ALL CJ*S FOR SINGLE XI 
IF(.NOT«CUVSE) GO TO 180 
c 
Bare ee HcRE FOR LAMBDA’S JARYING FIRST THEN XI*S 
c 


95 CONTINUE 
0O 160 LE=1,NE 
OO 150 LL=1i,yNL 
CULL) =H=4e SL AMBOA(LL?) 
Ou 110 LS&=1,NSE 
OCG 180 LSEA=14,NSE 
C(LSE,LSEA) =A(LSEyLSEA) tLAMBDAC(LL)*B(LSE,LSEA) 

100 CC(LSE,LSEAD=C(LSEsL SEA) 
E(tSe) =ePS(lLe LSE) ¢LAMBDA (LL) FFT (Le,ylLSe) 

110 SE(LSEV=E(LSE) 

C#sere CALL SIMULTANEOUS EQUATION SOLVER 

CALL SLEQ(CycyNSEy IER) 
IF CLeERscQeit)d WRITE (6,1008) 
SUMIE(LL) =0-6 
O00 120 LSE=1,NSE 
EACLE LSE) =EC(LSE) 
IFCFLAG) DOT(LLyLSE) =E (LSE) 
SUMIE(LL) =SUMIE(LLIFINILSEDFECLSED 
IFCFLAG) OS(LL)=SUMIECLL) 
IF(.eNOTeBETAC) GO TO 120 
PLCLA=VAR(LE) 


RA=1L eS (KAPPAF( (te /PLCLA) #¥(PLOLAS(PLELAW16) P= (1 eo /PLOLA) ** 


1(1e/(PLCLA-1.)))) 
SACLL yg LSED=(RAZ4.)*CA(LL LSE) 
120 CONTINUE 


IF(BETAC) SUMIE(LL)=-RA¥ ( (ie PLSLAD* IL}. 25*SUMIc (LL)? 


IFCOIVSR(LE).SQe1+) GO TO 140 
00 130 LSE=1,NSE 
130 VOP(LLE LSE) ec (LSE) /O0ITUSR(LE) 


_ BIST AVAILA 


BLE COPY 
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PROGRAM MAIN 


115 


120 


125 


130 


135 


140 


145 


150 


155 


160 


165 


170 


SOP (LL) =SUNTe (LL) ZOIVSRILED 
140 CONTINUE 
G#eeSE IF OPTION IS TRUE EXTeNOcCO OUTPUT IS GIVEN 
IFCOPTION) CALL WRITE (LeyLL) 
150 CONTINUE 
Crrere IF OPTION IS FALSE STANDARD OUTPUT IN FORM 
C¥Fere ALL CJ*S FOR SINGLE XI IS GIVEN 
IF(.NOTe«OPTION) CALL WRITS(LE,NL) 
160 CONTINUE 
IF(.NOTeFLAG) GO TO 280 
FLAG=.FALS=. 
NE=0 
CUdSt=eFALSee 
GO TO 10 
180 CONTINUE 


c 
ceeere HcRE FOR STANDARD OUTPUT IN FORM 
Ceeaee ALL XI*S FOR SINGLE CJ 
(m 
D0 260 LL=1,NL 
CUCLL) =46/LAMBOACLL) 
00 250 LE=1,NE 
00 210 LSE=1,NSE 
DO 200 LSEA=1,NSE 
C (LSE yLSEA) =A (LSE yLSEA) +LAMBOA(LL)*B(LSE,LSEA) 
200 CCLLSEsLSEA)=C(LSEyLSEA) 
2 (LSE) =cPS(LE,LSE) HL AMBDA (LL) *FT (LE, LSE) 
210 EB(LSE)=E (LSE) 
crests CALL SIMULTANEOUS EQUATION SOLVER 
CALL SLEQ(C,EyNSE, IER) 
IFCIERsEQe1) WRITE(6,1008) 
SUMIE(LE) =0. : 
00 220 LSE=1, NSE F 
EA(LEyLSE) =E(LSE) 
220 SUMIE(LE) =SUMIE(LE) tKIN(LSE)*E (LSE) 
IF(DIdSR(LEDeEQe16) GO TO 250 
DO 230 LSE=1,NSE 
230 DOP (LeyLSE) =E (LSE) ZOIVSR(LE) 
SOP (LE) =SUMTé (LE) /DIVSR(LE) 
250 CONTINUE 
IF(.NOTsCMPHAI) GO TO 258 


c 
Creer HERE FOR CALCULATION OF H(C)*S AND SCRIPT I*S 
c 

NEM4G=NE<4 

OG 254 LSE=1,NSE 

HC(LSE)=06 


00 253 JN=1,NEM4 

IF(JNeNE/1) GO TO 252 

YO1=Y (JN) *O00P (JN4+15LSE) 

YO2=Y (JN#1) *ODP(UNt2,L SE) 

YOZ=¥ (UN#2) *DOP(IN+3 ,LSE) 

X21=VARCINF2) -VARCING1) 

X SZ=VARCINt3) -VARCING2) 

XSL=VARCING3) -VARCING1) 
HCINTR=(X21/60) *(60* YO1+ (204 (X31/X32)) *CVO2-YD1) = (X21 * #2 
AZ (XZL*X32)) *LYOS-YO1)) 


, _ BEST AVAILABLE COPY 
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PROGRAM MAIN 


c 


253 


254 


255 


256 


GC TO 253 

YO1=¥ (JN-1) *Uu0P(JN,LSE) 

YO2=Y (JUN) *UOP (UN+1,LSE) 

YU3=Y (JN¢1) *DUP(UNt2,LSt) 

KX12=VAR( JUN) -VARCINGL) 

X1S=VARCIN) -VARCUN#2) 

X32= JARCINt 2) - JARCINt1) 

HCINTR=(X 32/60) *( 60% YO2+ (2+ 0 X12/KX13) )* (YO3-YD2) = (X32*%*2 
AZ (X12*xX13))*(Y01-¥02)) 

HC(LSEVSHCC(LSED FHCINTR 

KTLESCONST(1) *2A(1,LSE)+CONST (2) *EA(2,LSE) 
KHLE=CONST(3) *cA(1,L SE) +CONST (4) FEAL 2,LSE) 

HLe (i Sc) =KTLe #¥ (1) ¢KHLE* YPLE 

KTTL=CONST(5) *CA(Ne=25LSc) tCONST (6) *EA CNEHW1,LSE) CONST (7) *EACNE, 
ALSL) 

KHTE=CONST (8) *ZA(Ncm2,LSE) tCONST (9) FEA CNE@1,LSc) ¢CONST (10) *E ACNE, 
ALS2) 

HTE (LSE) SKTTE*¥ (NE M3) #KHTE*VPTC 

H(LSt) =HLE(LSE) +HC (LSE) ¢HTECLSE) 

CONTINUE 

ISPTC=6. 

DO 256 JN=1yNEM4 

IFCUNeNE21L1) GO TO 255 

YS1=Y (JN) *SOP (JNt1) 

YS2=¥ (UNt1) *SOP(CUNF2) 

YS3=¥ (CUN#2) *SOPCINF3) 

X21=VARCINt2) -VAR(IN#1) 

X32=VARCIN+3) -VARCIN#2) 

X31=JARCIN+3) -4AR(ING1) 

LINTR=(X 21760) (60 FYSL 4020+ 0 X31/X32) DECYS2*V¥S1) = OX21%*2 
AZ (X31"X32)) #CYS3-YS1)) 

GO TO 256 

YS1=V(JN=-1) *SOP(JUN) 

YS2=Y¥ (JUN) *SOP (JN+1) 

YSS=Y (UN#1) *SOP(UN#2) 

K12=VAR( UN) =JARCUN¢1) 

X13S=VAR( JN) “VARCIN#2) 

X32=VAR(IN+2) -VAR(CING1) 

LINTR=(X32/ 606) (60 FVS24(20¢(X127K13) DFC VS3=¥S2) = (X32**2 
AS OX12*X13)) FL YS1-YS2)) 

ISPTC=LSPTC#LINTR 
SKTLE=CONST (1) *SUMIE (1) #CONST (2) *SUMIE (2) 
SKHLE=CONST (3) *SUMIE (1) +CONST (4) *SUMIE (2) 
ISPTLE=SKTLE*Y(1) *SKHLE* YPLE 

SKTTE=CONST (5) *SUMIE (Ne -2)*#CONST (6)*SUMIc (NE=1)+CONST (7) *SUMIE (NED 
SKHTE=CUONST (8) *SUMIE (NE=2) #+CONST (9) *SUMIE (NE=1) + CONST (10) *SUMIE (NE 
1) 

ISPTTE=SKTT&E* Y(NE=3) tSKHTE*YPTE 
ISCRPET=ISPTLE+ISPTC+ISPTIT= 

IF(.NOTe«BETAQ) GO TO 256 


Gers te H=oRE FOR QUADRATURE METHOO BETA CALCULATIONS 


c 


YPILOS=YPT=*(IL+OS (LL) ) 
OCM=-YPILOS+#ISCRPT 
B1ISUM=0. 

O00 257 LSe=1,NSE 


BEST AVAILABLE COPY 
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PROGRAM MAIN 


YPUD (LSE) =Y¥PTE*OOT(LLyLSc) 
SETA(LSe) =YPDO(LSE)-H(LSE) 
BISUM=bISUM#BETA(LSE)*IN (LSE) 
OCMASK=-YPTE*IL-8ISUM 


CONTINUE 
CALL WRITECNE yLL) 
CONTINUc 
CONTINUE 


FORMAT (8H IcR = 1) 

FORMAT (1043) 

FORMAT (A5 55X9A10,A595X9N555X_A5) 
FORMAT (1X 9 A559 5X9 A109 A5595X 9 N555XyA5) 
FORMAT(1X,10A8) 

STOP 

ENO 


BEST AVAILABL: Co? 
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SUBROUTINE WRITE 


SUBROUTINE WRITEC(LEyLL) 

CGOMMON/ALL/A (959) 934999) 950999) pCO(999) E09) EBL 9) ,CU 9) 

1, EA (3059) EPS (3059) p FT (3059) »HEAD (30517) »DOP (30, 9), SUMIE( 30), S0P(3 
20) ,DIVSRI30), VAR(30) sCUVSi pL AMBIA(30) »NSEyNLyNey OPTION,FIRST 
3,HLC9) sISPTC, ¥ (30) pCMPHATI pHLE (3) gHTE (9) pH( 9), LSPTLE,ISETTE,ISCRPT 
4»hEADE(S) sPLCLA,HEADH( 10) »8c TAQ, dETAC, BETA(9) , YPDO(9) »FLAG,OCM 

Ds KaPPA,RA,OS(30) ’ OCMASK, YPILOS 

INTEGER HEADE yPLCLAy HEADH 

LOGICAL OPTIONsCUVSE,FIRST sCMFHAI sBETAU,BETAC,FLAG 

REAL LAMBDA,KAPPA 

ReAL ISPTC,ISPTLe,yISPTTE,ISCRPT 

c 
NEM3=NE<3 


cr 
v 


Cerere TEST FOR cXTENOED OUTPUT OPTION 
c 
IF (.NOT.OPTION) GO TO 200 
IFC(.NOTeFIRST) GO TO 130 
FIrST=eFALSE. 
TH=0 
c 
COeere PRINTS A(MyN)*S ANDB(M,N)*S (ONLY ONCE) 
c 
WRITE (o,1010) 
OC 110 LS&=1yNSE 
110 WAITE(6,1083) (ACLSE,LSEA) »LSEAH4,NSE) 
WRITS (6,1011) 
OC 120 LSE=1,NSE 
120 WRITE (56,1005) (B(LScyL SEA) » SEA=1yNSE)D 
c 
Ger tee HERe TO PRINT OUTPUT IN EXTENDED FORM 
c 
130 IFC2*(LL/2) NEoLL) WRITE (53,1095) 
WRITE (65,1102) 
WRITi (691191) (HEAD(LEyI) »1=1517) » 
ARIT2(6,1099) VAR(LE) yDIVSRILE) »LAMBDA(LL) ,CU(Lid 
WRITE (6,1098) (EPSCLEyLS=) pL SE=1y NSE) 
WRITE (6,1097) (FT(Le,LSE) »LSC=1,NSE) 
ARIT= (691090) (EB(LSE) »LS==1,NS=) 
WRIT=(o,1012) 
WRITE (60,1003) C(CO(LSE,LS=A) yLSEA=1y NSE) y LSEH1yNSED 
WkITE (691013) (HEAVE (1) »L=159) 9 (HCADE(T) yI=2,3) 
WRIT2 (691050) CEACLL,) LSE) »LSt=1yNSc) »SUMIE(LL) 
IF (COIVSR(L)efQei.) GO TO 159 
Wk1Tc (610106) 
WRITE (651030) (DOP(LLyLS=-) »LSE=1, NSE) »SOP(LL) 
CONTINUE 
RETURN 


CONTINUE 

TFC eNOT.FIRST) GO TO 204 
IF(FLAG) GJ TO 205 
FIRST=eFALSE. 

IF (eNOTeCMPHAI) GO TO 205 


atch HeRe TO PRINT ALRFOIL CAMBER LINE COORDINATES 


BEST AVAILABLE Copy 


Lui 


AFFDL-TR-77-63 


SUBROUTINE WRITE 


60 


65 


70 


75 


60 


65 


90 


95 


100 


105 


110 


203 
204 
c 


Crteee 


c 
205 


c 


cereee 


c 


207 
210 


220 
250 


Cc 
Cc 


Crane 


c 
300 


310 


340 


WRITE (6,1095) 

WRITE (6,1001) (HEAOH(I),1=1,10) 
OO 203 J=1,NEM3 

WRITE (6,1002) VAR(J+1) ,¥ (J) 
IF(BETAQ) GO TO 400 


HERE TO PRINT STANDARD OUTPUT 


CONT INUE 
WRITE (6,1095) 
IF(.NOTeCUdsSE) GO TO 300 


HERE FOR ALL CJ“S FOR SINGLE XI 


WRITE (6,1101) (HEAD(LE,I) »I=1,17) 
WRITE (651103) VARCLE) ,OIVSRILE) 
IF(BETAC) WRITE(6,1128) RA 

WRITE (6,1098) (EPS(LE,LSE) »L. SE=1, NSE) 
WRITE (6,1097) (FTC(LE,LSE) pLSE=1,NSE) 
IF(eNOTeBETAC) GO TO 207 

HEAUE(1)=5H BETA 

HEADE (2) =10H 
HeADE (3) =5H 
WRITE (651114) 
00 210 N=1,NL 
WRITE (651033) CUCND» CEACN LSE) gL Sc=1yNSE) »SUMIE(N) 
IFCOIVSR¢LEd+EQs1e) GO TO 250 

WRITE (6, 1035) 
WRITE (6,1101) 
WRITE (6, 1103) 
WRITE (651098) 
WRITE (6) 1097) 
WRITE (6,1112) 
00 220 N=1,NL 


OCM 


(HEADE (1) 12159), (HEADECT) »I=2,3) 


(HEAD(LE,I) »I=1517) 
VAR(LE) »OIVSRILE) 
(EPS(LE LSE) »LSE=i,NSE) 
CFT CLE LSE) »LSE=1, NSE) 


WRITE (651033) CUCN), (OOP (NyLSED pL SE=12yNSE) ySDOP (ND 
CONTINUE 
RETURN 
HcRE FOR ALL XI*S FOR SINGLE CU 
WRITE (691101) (HEAD(LE I) y»I=1517) 
WRITE (651107) CULL) 
WRITE (6,1115) (HEADE (1) » 13159), (HEADECI) y 12,3) 
00 310 N=1,NE 
WRITE (651033) JARCN) » CEACN,LSE) »LSE=1) NSE) »SUMIE (ND 


IF(OIVSR(LE) cEQeie) GO TO 380 

WRITE (6, 1095) 

WRITE (6,1101) (HEAD(LE,I) »1=1,17) 

TFC HEADE (1) eEQe5H SETA) WRITE (6,1125) 

IFQ(CMPHAL) WRITE (06,1004) 

WRITE (6,1107) CJCLLd 

WRITE (6, 1113) 

00 340 N=1,Ne 

WRITE (651033) VAR(N) » (OCOP(N, LSE) yLSE=1yNSE) ySOPCN)D 
TFCeNOTeCMPHAI) GO TO 389 


. BEST AVAILABLE 


ran 


Wes i 


tt cee er etree at 
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SUBROUTINE WRITc 


115 


120 


125 


130 


135 


140 


145 


150 


155 


169 


165 


170 


5 pt A ROR 


c 


Crreee 


Cc 


380 


400 


410 


500 


WRITE (6,1120) 
WAITe (6,1121) 
WRITt (0,1122) 
WRITE (6,1123) 
WRITE (60,1124) 
CONTINUE 
RETURN 


CONTINUc 


CHLE (LSE) »LSE=1,NSC),ISPTLE 
(HC (LSE) pL Sc=1_NSc) pISPTC 
(HTE (LSE) »LSE=4y,NSE),LISPTTE 
(H(LSE) »LSE=1) NSE), ISCRPT 


HeRE IF PRINTING ONLY HyI1SCRPT, YPOO,BETA,ANO DCM 


IF CCIH/3) *3-NEeTHeANDeIHeNEo0) GO TO 410 


WRITE (0,1095) 
WRITE (06,1129) 
WRITE (0, 1004) 
WRITE (o, 1000) 
CONTINUc 
IH=IH+1 
WRITc (0, 1117) 
WRITE (6,1120) 
ARITE (6,1121) 
WRITC (60,1122) 
WR1Tc(o, 1123) 
dk1Te (o,1124) 
WRITE (6,1126) 
WKITct (60,1127) 
WAITe (6,1130) 
CCNTINUE 
RETURN 


(HEADH (I), T=1,10) 


CULL) 


CHLE (LSet) yLSES1,NSE) »LSPTLE 
(HC (LSG) yLS=E=1yNSE) »ISPTC 
(HTE (LSE) yLSES1,NSE),ISPTTE 
(H(LS=) »LSe=1yNSe) y ISCRPT 
(YPOO(LSE) yLSc=1yNSE) y YPILOS 
(BETALLSE) »LS==1,NSE) DOM 


OCMASK 


999 FORMAT (1X 510A38) 


1060 


FORMAT (1043) 


1001 FORMAT(26X,10A6/ 
A50X,S1HAIRFOIL CAMBER LINE COOROINATES/ 
B52X,27HFOR HC RANGE OF INTEGRATION// 
CO7TXyLHXy1OXy1HY) 


1062 
1003 
1004 
1010 
1011 
1012 
1013 


FORMAT (48X9F 14099 5XyF 1207) 


FORMAT (1Xy9F 1206) 


FORMAT(/26H HC RANGE XI= .05 To 


0975) 


FORMAT (7H1A (999) 9 104Xy1A1094X,1A10,7) 
FORMAT (//7H B(999) 7/) 
FORMAT (23H ACI,J) + LAMBUA*B(I,J)7) 

FORMAT (/4X 9A5 9 3H(0) yp SX A5y SHULD yp SXyASy SHV 2) pSXyASy SHS), 


ASX 9 Ady 3H (4) 95X9A5y 3H (5) 9 5X 9A 5_3H( 6) y SX 9A5y 3HI7) 9 SX9ASy SHU8) 


B4X,A10,A5/) 


1016 FURMAT (/6X,6HDOP(0) » 7X yGHOOP (1) 5 7Xy6HOOP (2) 97Xy5HOOP(3),y 


1030 
1033 
1095 
1096 
1097 
1098 
1099 


A7Xy6HOUP (4) 57X,6HOOP (5) 9 7X y6HOOP (0) » 7X y6HOOP (7) 5 7XyG6HIOP( A), 


B12XySHSOP/) 


FORMAT (1X9 9F1 307, 3X, F L307) 


FORMAT (1X 9F 90 592Xy9F 1207 5 4XyF 1166) 
FORMAT(LH1,110X91A1054Xs1A109/////) 
FORMAT(/19H ET(M) +LAMBDA*FT(M) 9/1X,9F 12657) 
FORMAT(/6H FT(M) s9F1265) 
FORMAT (/6H ET (M)y9F12.5) 


FORMAT (1% s14H (XI) OR 


(A) 


63 


SyF120e7» 


BEST AVAILAB 


IZ CAD 


iLL 


LUE 


if 
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SUBROUTINE WRITe 


& 28Xy7HDIVSR =9F120¢395X910H LAMBDA = 9F12.5,10X96H CU = yg 
AF 12.5) 

1101 FORMAT (1X,17A8) 

L102 FORMAT (AX yg L3IOHEFE te HEHEHE EEE HET EE EE EEE HFEF EEE HEH Ete et ett eee ee tet et 
AT EEE EET EE EE EE EE EE EEE EE EET EEE TEE EE EEE EE EEE EEEE EEE SHEESH ETE EET EE Hee tt 
BDeteee+eeeeeeeetetse) 

1103 FORMAT(10X,9HXI OR A =yF12069510X,8HDIVSR = »F12-.9) 

1107 FORMAT(SX,5HCU = ,F10.€) 

1112 FORMAT (/0X52HCJ 99X »GOHDUP (0) y OX, OHOUP (1) » 6X» 6HOOP (2) 6X, 6HODP( 3) 
AoX,OHODP (4) 6X,6HDOP (5) y OX» 6HDOP( &) » OX »OHOOP(7) » 6X,» 6HO0P( 8) 4 
810X,4H SOP/) 

1113 FORMAT (/6X,2HX1,9X,6HOGP (0) »6Xy6HODP (1) »6Xy6HOOP (2) »OX,6HOOP (3), 
A6X,6HDOP (4) »6X,6HDDP(5) 6X y6HODP (6) 9 6X »OHDOP (7) »6X,6HDUP (8), 
B10Xs4H SOP/) 

1114 FORMAT (0X s2HCSy7XyA5y3H(0) yy 4XyA5y SHO1) 9 YX gA5y SHU 2) 9 4X yATy SHIT) 
Ny 4X y AD 3H 14) 9 4X 9A5y3H(5) p4XyA5y 3H (6) p4X ep AS yg SHIT) 9 4XyA5 5 3H(8) y 
B3X,A10,A5/) 

1115 FORMAT (/6X_ 2HX1 9 7X 9A5y 3H(0) py 4X gAS yg SHI1) p4XyASy SHW2) 9 4X yASy SHO 5S) 
Ay GX» AS,3H(4) 9 4X,A5 93H(5) 9*XyAS 93H(6) 9 4X 9A5e3H(7) 94X,AS 93H(8) ’ 
B3xyA10,A5/) 

1117 FORMAT(/1X5HCU = 4 F106 6) 

1120 FORMATI/10X ,4HHOO) 58X y4HH(1) y 8X9 YHHO2) py BX yp 4HH(3) 5 BX, HHHI a), 
ABXs%HH(5) 93Xy 4HH(O) 9 BX pGHH(7) » BXy HHH (8) y 14Xy6HISCRPT) 

1121 FORMAT(GH Lky9F1207,3X,3HLE yF13.6) 

1122 FORMAT(GH C 59F120753Xy3HC F136) 

1123 FORMAT(GH TE 9 9F120753Xy SHTE 9F 1326) 

1124 FORMAT(GH TOT ,9F126753Xy3HTOT F135 0657) 

1125 FURMAT(1X,9HKAPPA = 1) 

1126 FORMAT (9X »7HYPUD(D) SX e7HYPDD(1) »5X»y 7HYPDD(2) 5X y7HYPDO(3), 
A5%y7HYPDD (4) 9 5X9 7HYPOD(5) 9 5X9 7HYPID(6) y5Xy7HYPDD (7) 95Xy 
B7HYPDO (8) 51 2X yGHYPIL OS s/4Xy9F 120 7,6X9F1366/) 

1127 FORMAT(9X y7HBETA(0) »5Xe7HSETA(1) 9 5X» 7HBETA(2) 5X y7HBETA(3),y 
ASX y7HBETA (4) 5 SXy7HBETA(S) 5X y 7HBE TACO) 5X» 7HBETA (7) 9S Xy 
B7HEETACS) 5 14X y SHOCM/ 4X9 9F 120 7y6Xy F136) 

1128 FORMAT(/10X,4HRA =,F 12-7) 

1129 FORMAT(39H H(N) ANDO Bc TA(N) COEFFECIENTS AND SUMS ) 

1130 FORMAT (116X,F13.6,1H*) 

END 
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(a) REGULARLY BLOWN FLAT PLATE 
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(b) SINGULARLY BLOWN FLAT PLATE 
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(c) REGULARLY BLOWN,MECHANICALLY FLAPPED AIRFOIL 


FIG 1. Fundamental Jet-Flapped Airfoil Cases 
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Vertical Scale is Exaggerated 


Jet Sheet 


FIG 2. Cambered Jet-Flapped Airfoil 
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Superposition Principle 
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FIG 4. Notation for Quadrature Formula 
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FIG 5. Normalized Power-Law Camber Lines 
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FIG 6. Interference Lift and Pitching-Moment Coefficients 
for Basic Power-Law Camber Lines of a = 2, 4, and 
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